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ABSTRACT 

Context. Anti-proton and positron Galactic cosmic ray spectra are among the key targets for indirect detection of dark 
matter. The boost factors, corresponding to an enhancement of the signal, and linked to the dumpiness properties of the 
dark matter distribution, have been taken as high as thousands in the past. The dramatic impact of these boost factors 
for indirect detection of antiparticles, for instance with the PAMELA satellite or the coming AMS-02 experiment, asks 
for their detailed calculation. 

Aims. We take into account the state-of-the-art results of high resolution N-body dark matter simulations to calculate 
the most likely energy dependent boost factors, which are linked to the cosmic ray propagation properties, for anti- 
protons and positrons. The results from extreme, but still possible, configurations of the clumpy dark matter component 
are also discussed. 

Methods. Starting from the mass and space distributions of sub-halos, the anti-proton and positron propagators are 
used to calculate the mean value and the variance of the boost factor for the primary fluxes. We take advantage of the 
statistical method introduced in Lavalle et al. (2007) and cross-check the results with Monte Carlo computations. 
Results. By spanning some extreme configurations of sub-halo and propagation properties, we find that the average 
contribution of the clumps is negligible compared to that of the smooth dark matter component. Dark matter clumps 
do not lead to enhancement of the signals, unless they are taken with some extreme (unexpected) properties. This result 
is independent of the nature of the self-annihilating dark matter candidate considered, and provides precise estimates 
of the theoretical and the statistical uncertainties of the antimatter flux from sub-halos. 

Conclusions. Spectral distortions can still be expected in antimatter flux measurements, but scenarios invoking large and 
even mild dumpiness boost factors are strongly disfavoured by our analysis. Some very extreme configurations could 
still lead to large enhancements, e.g. (i) very small clumps with masses < 10~ 6 Mq following a M~ a mass distribution 
with a > 2, highly concentrated with internal r~^ profiles with f3 > 1.5, and spatially distributed according to the 
smooth component; or (ii) a big sub-halo of mass > 10 7 A/q within a distance of < 1 kpc from the Earth. However, 
they are very unlikely from either theoretical or statistical arguments. 

Key words. Dark Matter 
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1. Introduction 

The existence of dark matter (DM) has been established 
by various astronomical observations, from galactic to cos- 
mological scales. The evidence come from gravitational ef- 
fects, such as the observation of the rotation curves in spiral 
galaxies and velocity dispersion in elliptical galaxies, the X- 
ray emission and peculiar velocity dispersion of galaxies in 



Send offprint requests to: 
Iavalle@in2p3 . f r 
yuanqOmail . ihep .ac.cn 
dmaurin@lpnhe . in2p3 . f r 
bix j Omail . ihep .ac.cn 



the clusters of galaxies and the weak lensing effects, all in- 
dicating much steeper gravitational potentials than those 
inferred from the luminous matter. Recently, there have 
been two stro ng smoking guns f rom the Bullet clust er sys- 
tem 1E0561 ([Clowe et alj|2004 120061 iBradac et ai1l2006h 
and a DM ring discovered around the cluster CL0024+17 
(|Jee et al.ll2007| ) , which may indicate the existence of DM in 
the sense that it has first provided means to study the dy- 
namics of DM itself. Note, however, that mo dified gravity 
models might still offer a viable alternative (lAngus et al.l 
120061 . 12007t lAngus fc McGaugh l l2007t iFamaev et al.ll2007D . 



The nature of DM is still unknown, remains one of the 
most outstanding puzzles in astrophysics and cosmology, 
and is challenging from the particle physics view point. 
Nevertheless, the unprecedented precision reached in obser- 
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vational cosmology in the last decade, thanks to the com- 
bined use of different probes (CMB, type 1A supernovae, 
large scale structures, deep surveys, primordial abundances, 
etc.), yields a rather precise estimate of the total amount of 
non-relativistic matter in the Universe, encompassing the 
standard baryonic matter, of which density can be pre- 
dicted and measu r ed independently ( for reviews, see e.g . 
iTvtler et all |2000| ; iLesgoureued 12004 lYao fc et all 120061) . 
The overall contribution of matter to the (critical) energy 
density of the Universe is ~ 30%, while the baryonic compo- 
nent accounts for 4% only. Hence most of the matter should 
be dark and of non-baryonic origin, requiring physics be- 
yond the standard model of particle physics. The most at- 
tractive scenario involves weakly interacting massive parti- 
cles (WIMPs). An appealing idea is that WIMPs could be 
thermal relics of the early Universe, which naturally give 
rise to a cosmological abundance in the range of the ob- 
served value if both the interaction strength and the masses 
are taken at the weak scale. Indeed, because of their ther- 
mal origin, WIMPs should still (weakly) interact with ordi- 
nary matter, and even annihilate if they are preserved from 
matter-antimatter asymmetry. Such particles can originate 
naturally in the context of supersymmetric (SUSY) or 
extra-dimensional (ED) extensions of the standard model, 
independently developed to tackle the issues of the unifica- 
tion of interactions and energy scale hierarchical problems. 
Indeed, in such theories, the stability of the proton is very 
often ensured by the conservation of some new discrete sym- 
metry that guaranties the lightest exotic particle to be sta- 
ble. Such paradigms provide very good candidates for DM 
(for reviews, see e.g. lJungman et all 1996; Bergstrom 20001 : 
iBertone et al.ll2005al) . In particular, the minimal SUSY ex- 
tension of the standard model (MSSM) can yield DM par- 
ticles, the most famous being the neutralino, a Majorana 
fcrmion. The cosmological constraints on the SUSY pa- 
rame ter space have been extensively studied in the litera- 
ture dBelanger et all 120061 : iDiouadi et alll2006l : iBaer et all 
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present running or proposed experiments, either directly by 
measuring the recoil ene rgy when they scatter off a detec- 
tor nuclei ((Munoz 2004), or indirectly by observing their 
annihilation prod ucts, such as anti-pro t ons, positrons, 7 - 
rays or neutrinos (|Bertone et al.l [2005at ICarr et aTll2006h . 
They may also be generated in the next generation collid- 
ers, which is the most direct way to probe the existence of 
new particles. The direct and indirect detection methods 
are viable and complementary to collider studies in order 
to further constrain the nature of DM. 

For indirect detection in the Milky Way, since the anni- 
hilation rate is proportional to the square of the DM den- 
sity, the Galactic Ce ntre is believed to be a p romising source 
of DM annihilation (jBergstrom et al.lll998| ). However, the 
existence of the central super-massive black hole and the su- 
pernova remnant Sgr A* are likely to heavily contaminate 
the DM s ignals with high-energ y standard astrophysical 
processes (lAharonian et al.ll2006D . Alternative sites, such 
as the DM dominated dwarf spheroidal galaxies (dSph) 
orbiting close around the Milky Way, or even DM sub- 
structures inside the Milky Way, could be more favourable. 
Indeed, the existence of a myriad of sub-halos throughout 
galactic-scale host halos is a generic prediction of the cold 
dark matter (CDM) paradigm of structure f ormation in the 
Universe. H igh resolution simulations (e.g. iDiemand et alj 
show that for the ACDM scenario, the large 
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scale structures form hierarchically by continuous merg- 
ing of smaller halos. As remnants of the merging process, 
about 10% to 50% of the total mass of the halo could 
be in the form of sub-halos. Moreover, the centres of sub- 
halos, like their hosts, are found to have high mass densities 
and therefore, could be ideal targ ets for 7-rays searches of 
W IMP annihilation prod ucts (e.g. lKoushiappas et al.|[2004l : 
IbI 120061 iBi et al.l [20071 and references therein). A long- 
standing issue is the possible overall enhancement — boost 
factor — of the signals from the smooth component, due 
to th e presence of such inhomogeneities (jSilk fc Stebbinsl 
Il993l ). The first studies dedicated to indirect detection of 
DM focused essentially on 7-rays, and more marginally on 
anti-protons, but suffered fr om the lack of informa tion on 
DM sub structures (see e. g. Bergstrom et al.|[l999l ). More 
recently, iBerezinskv et al.l (|2003l ) discussed in more details 
the 7-rays case, finding boost f actors no larger t han a few. 
Furthermore, a recent study bv lDiemand et all (|2005af ) re- 
heated the debate on dumpiness, because the authors, by 
means of a very high resolution N-body experiment (but 
stopping at z = 26), found that the Galaxy could be popu- 
lated by a huge number density of sub-halos as light as the 
Earth. While the survival of such light clumps against tidal 
effects is still questionable, they could yield a significant 
contribution to the Galactic diffuse 7-ray flu x by assuming 
a very cuspy sub-halo profile (|Bi et a l. 2006). Nevertheless, 
some recent works also indicate that the current parame- 
ter ran ge for dumpiness may provide only marg inal global 
effects (|Diemand eraT1l2007al iPieri et al.ll2007l) . The aim 
of the present paper is to provide a detailed study of the 
impact of cosmological sub-halos on the primary antimat- 
ter Galactic cosmic ray (GCR) flux, as elaborate as that 
already performed for 7-rays. 

In lMaurin fc Tailletl (|2003f ). the authors noted that the 
difference in propagation properties for p and e + was likely 
to translate i nto different boost f actors for these species. 
More recently, lLavalle et al.l (|2007f ) provided a detailed for- 
malism to tackle the calculation of antimatter CR fluxes, 
when boosted by DM dumpiness. They showed how the un- 
certainty on the spatial distribution of clumps transfers to 
an uncertainty to the predicted boosted cosmic ray positron 
flux, an effect that depends on energy. More generally, this 
effect depends on the clump number density in a volume 
bounded by the characteristic diffusion length of the in- 
volved species. For the sake of clarity, these authors have 
used a very simple model, in which all clumps have the same 
internal properties (masses and intrinsic luminosities), and 
mainly stressed the effects coming from their space distri- 
bution. Using this method. iBrun et al.l ((2007) fully treated 
a particular class of DM inhomogeneities— the intermedi- 
ate mass black holes (jBertone et al.l [2005bh — finding large 
boosts with huge variances for the signals: such large vari- 
ances tag unpredictive scenarios. This means that in the 
case of a positive detection, such scenarios can certainly be 
tuned to reproduce the data, but generally at the cost of 
a vanishingly small associated likelihood for this configura- 
tion. 

I n this paper, we stud y a more natural DM scenario 
(e.g. IDiemand et a l. 2005a), in which substructures fill the 
whole Galaxy down to a minimal mass M m [ n > 10~ 6 M Q , 
with a mass distribution dN/d\ogM oc M~ Um (a m ~ 0.9), 
and a cored spatial distribution. We survey different DM 
configurations in great details by using different sub-halo 
inner profiles, different mass distributions or different con- 
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ccntration models (this has alread y been well stud ied in 
the context of gamma-rays, see e.g. lUllio et aill2002[ ). It is 
important to better quantify the boost and v ariance of anti- 
matt e r signals since the sa tellite PAMELA (jPicozza et al.l 
l2006t ICasolino et al.l [20071 ) . successfully launched in June 
2006, will soon provide new results on antimatter fluxes. 
The DM description suffers uncertainties, and its impact 
on the calculated fluxes adds up to th e existing uncertain- 
ties fr om the propagation p arameters (jDonato et al.ll200H 
12004 iLionetto et all 12005( 1. Regarding this latter issue. 
PAMELA should also update our current knowledge of the 
particles transport in the Galaxy, thanks to secondary-to- 
primary ratio measurements (e.g. B/C). This is crucial for 
the background calculation (standard antimatter produc- 
tion) in order to confirm/support any claim of an excess. 
Besides, AMS-02 should be launched in the coming years, 
and provide additional crucial information on GCR prop- 
agati on by measuring the radioactive species l|Battistonl 
l200l . 

Below, we take advantage of simplified for mulations for 
the p (e.g. iMaurin et al.ll2006a !) and e + (e.g. lLavalle et alj 
120071) propagators. Using the information of the mass and 
space distributions of sub-halos from N-body numerical 
simulations (see e.g . the r ecent Via Lactea simulation, 
iDiemand et all 12009 . l2007allbT ). we calculate the boost and 
the variance of the fluxes. We find that for all plausible 
choices of the clump properties and propagation parame- 
ters, boost factors for anti-protons and positrons are close 
to unity, with small systematic and statistical uncertainties. 

The paper is organised as follows. All relevant aspects 
(for this study) of the DM distributions in the Galaxy, in- 
cluding N-body simulation results are discussed in Sect. [5] 
The configurations retained are given in Sect.[3l where the 
key parameters entering the calculation of the clumpy flux 
(and its variance) are underlined. The propagation aspects 
are treated in Sect. [4j The methodology to calculate the 
antimatter flux, its variance and the corresponding boost 
factors is given in Sect. either by means of a semi- 
analytical approach (Sect. I5.2[) or by Monte Carlo (MC) 
simulations (Sect. I5.3|) . The reader not interested in the 
technical details can directly jump to Sect. where the 
results for positrons (Sect. 163)) and anti-protons (Sect. lO)) 
arc presented, highlighting the physical effects coming from 
clump properties, space distribution, mass distribution and 
GCR propagation. Because of the complex origins and the 
mixing of the relevant physical quantities, such details re- 
ally help to fully understand what kind of information boost 
factors actually encodes. We summarise and conclude in 
Sect. EJ 

2. DM distribution 

In the last few years, the advent of high resolution N-body 
simulations have increased the number of studies in this 
field, allowing for a better understanding and description 
of the DM dynamics. Even if many issues remain unclear, 
when comparing simulation results to the current obser- 
vations, collisionless codes now agree at the 10% level over 
wide dy namic ranges, providin g a robust framework for DM 
studies (jHeitmann et al.|[2007D . 

Throughout this paper, we will separate the WIMP an- 
nihilation contribution associated with sub-halos from that 
associated with a smooth component. The former will be re- 
lated to any DM inhomogeneity in the Galactic halo, inde- 



pendently of its physical scale — resolved or not in N-body 
simulations — while the latter will refer to the Galactic 
DM host halo itself which will be considered as a continu- 
ous fluid (again independently of the current resolution of 
N-body simulations). Although the Vlasov (or fluid) limit is 
likely to be reached when the number of particles involved 
in N-body experiments is huge, one should still be aware 
that such a statement is not trivial at all when dealing 
with the cosmological evolution of structu res, and that dis- 
crete ness might induce important biases (|Jovce &: Marcos! 
120071 ). Furthermore, one should also keep in mind that our 
DM modelling will rest on (or be extrapolated from) N- 
body experiment results, in the most precise of which the 
test particle mass is not lighter than ~ 10 4 Mq (cf. the Via 
Lactea simulation. IDiemand et a l. 2006), and for which the 
Vlasov limit is not reached at small scales. Nevertheless, we 
will assume throughout this study that the host halo pro- 
files of Milky- Way-like galaxies provided by N-body simula- 
tions describe a smooth fluid (WIMP gas), on top of which 
some sub-halos may be wandering. 

In the following subsections, we summarise the recent 
results concerning (i) generalities about DM distribution in 
halos of galaxies (Sect. l2TTj) and (ii) some specific consider- 
ations about sub-halo description (Sect. I2"?2"l and l2.3[ ). 

Given the scope of this work, we will merely consider 
spherical profiles. For sub-halos, several cases will be chosen 
to encompass some extreme (but still plausible) scenarios. 
This aims at providing realistic estimates of the boost factor 
uncertainties related to the clumpy DM component. 

2.1. Shape and profiles 
2.1.1. Spherical profiles 

A scale-invariant DM distribution based on N-body numer- 
i cal simulat ion results can be written in a general form as 
(|Zhaolll996D 



H {r/r s )->[l + {r/r s ) a ]^-'<)/ a ' V ; 

where p s and r s are respectively a scale density and a scale 
radius, which can be determined by measuring the relation 
between the mass of the dark halo and the concentration 
parameter from simulations. Such an empiric law can be 
used for galaxy cluster halos, galaxy halos and for sub- 
halos. In the following, wc focus on the central logarithmic 
slope a of the smooth halo component. Wc will discuss the 
Galactic scale radius r s and density p s in the section dealing 
with the concept of concentrat ion (cf. Sect. |2~. 2. 21) . 

Navarro, Frenk and White ([Navarro et al.lll997fl worked 
out the following set of parameters (a,/?, 7) = (1,3,1), 
which define the NFW profile, with a cusp scaling like r _1 
at radii smaller than r s . iMoore et alj (|l998l ) found another 
set with (o:,/3, 7) = (1.5,3,1.5) to fit their simulation re- 
sults, which is steeper than NFW at small radii, scaling 
like r -1 - 5 . More recent high resolution N-body simulations 
found that a NFW profile seems to underestimate the DM 
density in the central regions, while a Moore profil J] prob- 
ably overestimates it ([Navarro et al.l 12004 IDiemand et al.l 



Though those authors have since improved their early time 
results and parametrisations, we will still use these generic pro- 
file names to deal with inner shapes of profiles, for the sake of 
simplicity. 
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l2004al [2005bl and references therein). The mean slope of the 
cusp obtained from various co des is well fitted by a (1 , 3, 7) 
profile, with 7 = 1.16 ± 0.14 (jDiemand et al.ll2004aD . still 
in agreement with (7 = 1.3) analytical similarity solutions 
(|Ahn fc Shapirol 120051 ) . However, profiles may not have a 
universal shape (e.g. iNavarro et "all 120041 ; IStoehrl 12006ft . 
First, from the observational point of view, the relative scat- 
ter observed for the slope for fou r nearby low-mass spiral 
galaxies is 0.44 flSimon et alj|2005l ), three times larger than 
in simulations. Second, it was also recentl y stressed that 
asymptotic slopes may not be reached at all (INavarro et alj 



asymptotic slopes may not be rea cnea at all uiNavarro et ai.l 
2004lStoehrl)2006l; lMerri tt et al.ll2006l:lGraham et al.ll2006t 



Ricotti et alj|2007f ): according to iGraham et al.l ( 2006 ) . the 



Einasto function describes a simulated DM halo better than 
a NFW-like model. 

Closer to the Galactic centre, the super-massive black 
hole dominates the mass (r < tbh ~ 2 pc). The adia- 
batic growth of the black hole, if taking place in the cen- 
tre of the DM gravitational potential and without any 
merger, could lead to an enhanced DM density in this region 
(slope as steep as ~ 2.3 — 2.4, dubbed spike). Nonetheless, 
recent works seem to prefer a final r -15 behaviour for 
the DM density in the inner regions (see iMerrittl 12004 
iGnedin fc Primack1l2004l and references therein) . 

Finally, the luminosity of cuspy or spiky halos is sin- 
gular at the centre of the halo. However, a cut-off radius 
r cut naturally appears, within which the DM density sat- 
urates due to the balance between the annihilation rate 
\{<rv) p(r cut ) /w x ] _1 and the gravitational infalling rate of 
DM particles (Gp) _1/2 ijBerezinskv et alJll992h . Taking p 
about 200 times the critical density, we get 



2.1.2. Other open questions 



10- 



Mq kpc 



(2) 



During their history, structures undergo several merg- 
ers. The survival of the inner cusp of DM in these 
events has been investigated. The inner profile was 
found to be exceptionally robust, despite the relaxation 
that follows merging processe s (Boylan-Kolchin fc Ma! 
20041: lAceves fc Velazquezl 120061: IKazantzidis et al.l 120061 



McMillan et al.ll2007l : IValluri et al.ll2007ft . The implications 



are deep: the characteristic universal shape of the DM 
den sity profile may be set early in the evolution of ha- 
los (|Kazantzidis et aT1l2006l ). However, it is still not clear 
whether the central cusp is steepened or flattened when 
the baryonic distribution is taken i nto account. U s ing N - 
body hydrodynamical simu lations, Gnedin et al] (|2004ft . 
IGustafsson et alj (|2006f ) and lMaccio et all (|2006t ) find that 
the effect of gas c ooling steepens the inner density profile 
to 1.9 ± 0.2, while iMashchenko et all (|2006t ) claim that the 
random bulk motion of gas in small primordial galaxies 
(driven by supernovae explosions) removes the cusp, leav- 
ing only cored profiles for both small and large galaxies in 
the present Universe. 

Several other controversial issues remain and we only 
briefly quote them. The first one is the question of 
the halo evolution in the presence of a rotating stel- 
lar bar, leading to either a destruction of the cusp 
(see lSellwoodll2006l iMcMillan fc Dehnenl 120051 and refer- 
ences therein) or a steepening of the cusp ([Gnedin et al.l 



120041 : IColm et al.l 120061 ) . Some recent simulations includ- 
ing a stellar bar also hint at the emergence of a bar- 
like structure for the DM (DM bar) i n the central re- 
gion in the case of a strong stellar bar (| Colin et al.l 120061 : 
I Athanassoulal l2007h : this is the second issue, namely de- 
parture fro m sphericity. Direct observations eit her favour 
prola tness (lHelmill2004allbl: iLibeskind et~aT1l2005h or oblate- 
ness ([Johnston et al.ll2005t iMetz et alj|2007t iRuzicka et al.l 
I2007D . whereas for pure collisionless simulations , pro- 
latness is generally preferred (see also ITTett et al] l2007t 
iRomano-Diaz et al']|2007t iKuhlen et al.ll2007H 2 l . Prolatness 
for sub-halos is likely to depend on the position in the 
;alaxy, halos being more spherical in the outer regions 
Havashi et al]|2007l ). Then, more generally, there is some 
evidence that halos become more sp herical when the 
baryon i c cooling is taken in to account (IKazantzidis et al.l 
2004al: iNovak et al l 120061: IBerentzen fc Shlosmanl 120061 



Maccio et al]|2007btlDebattista et al.ll2007ft. or wh en a stel- 
lar bar is taken i nto account (IBerentzen "et~atll200l . or even 
during mergers (|Novak et al.ll2006ft . 



2.1.3. Simplifying assumptions 

It was shown that the choice of one or another DM pro- 
file for the smoot h component (jMaurin fc Tailletl 120031 : 
iDonato et al] 2004) is not crucial for the calculated flux of 
anti-protons and positrons. Indeed, charged particles dif- 
fuse on magnetic inhomogencities and fluxes are heavily 
suppressed (escape from the Galaxy) when originating far 
away from us, i.e. those from the Galactic centre. For cuspy 
profiles, the maximal difference is obtained between cored 
isothermal and Mo ore profiles (a factor < 2, see e.g. Fig. 2 of 
iBarrau et al.l l2005h . the difference between isothermal and 
NFW profiles bein g even smaller (< 20%, see Table II in 
IDonato et alj|2l304ft . Due to the lack of a definitive answer 
for the DM profile in the Galaxy (see the above-discussion), 
we will restrict ourselves in this paper to a spherically 
symmetric NFW profile (a,f3,j) — (1,3,1) for the galac- 
tic smooth distributions. Using triaxial halos or different 
profiles (e.g. 7 = 1.2 or any other profile) is expected to 
leave the main conclusions of the paper concerning the ef- 
fects of clump granularity in the Galactic halo unchanged. 
Sphericity is also assumed for the substructure^. Departure 
from spherical symmetry is left to a forthcoming study. 

Finally, we stress that although the existence of a DM 
spike in the Galactic centre is c rucial in the context of 7 - 
ray/neutrino indirect detection (jBertone fc Merrittll2005f ). 
its effect is merely not relevant in this study. This is due to 
the depletion of the signal through the diffusive transport 
of antiparticles, and also to the fact that GCRs originating 
from annihilations in the very tiny extent of this region are 
only a small fraction of the total yield that can reach the 
Earth (the DM annihilation contribution to the GCR flux 
is integrated over a diffusion volume instead of a line of 
sight for 7-rays). 



For further developments on the topic of tria xiality, asym- 
metrie s , as well as on the s pin of halos, see e .g . iMoore et al 
(|2004h : IZentner et alj (t2005h:]Lee fc Kand ll2006l):lGao fc White 
(200q); ICapuzzo-Dolcetta et all (|2006T ); IGustafsson et al 
(200$; iBett et all ([20071 1 . For the dependence of halo pa - 



rameters on the e nv ironment, see [M aulbctsc h et alj (|2007f ) ; 
lHahn et al.1 (|2007al fbl): iRagone-Figueroa fc Plionisl (|2007l ). 

3 Small struc tures, which formed earlier, are expect ed to be 
more spherical ([Moore et al.ll2004l ; lAlTgood et alJl2006f ). 
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2.2. Concentration parameter and sub-halo description 

The concentration parameter is a crucial quantity for com- 
puting the annihilation rates in (sub)structures. In this sec- 
tion, we actually present all relevant parameters that define 
a sub-halo. We will come back to the concentration (and the 
scale radius) associated with the host smooth halo just at 
the end (Sect. 

In the ACDM cosmology, the structures form hierar- 
chically bottom-up via gravitational amplification of initial 
density fluctuations. The properties of the emerging struc- 
tures and their subsequent evolution may be described by 
using the virial quantities. Foll owing the approach and def- 
initions of lBullock et al.l (|2001h . the two parameters p s and 
r Sl defined in Eq. ([T]), of a structure of mass M v - lr are ex- 
pressed in terms of the concentration c v i r and the virial ra- 
dius i? v ir- This outer radius is defined as the radius within 
which the mean density is A v ; r (z) times the matter density 
£l m Pc at redshift z. At z = 0, 



(47r/3)A vir (0)fW £ 



1/3 



(3) 



In the following, we use the standard ACDM values (fl m = 
0.24 p c = 148M(T) kpc" 3 ) and A vi r(0) « 340 ([Bullock et all 
120011: IColafrancesco et atll2006f ). 

The concentra tion parameter is defined 
(|Bullock et al.ll2001l) 



as 



Rvir 

r-2 



(4) 



where r_2 is the radius at which d/ dr (r 2 p(r)) \ r= r _ 2 = 
0. It was found ([Navarro et all Il997l : iBullock et ail l200lt 
lEke et al.ll200lh that c vir strongly correlates with M vlr , with 
larger concentrations observed in the first structures, i.e. in 
the lighter halos, which have formed in a denser Universe. 
This relation allows to express p s and r s in terms of the 
sole quantity M vlr . Note that the c- 



given for a NFW profile, for which 



vir 
fw 



M, 



rf w (M vir ) 



-Rvir (-Mvir) 
Cvir (Mvir) 



i r is generally 
2, so that: 

(5) 



This is easily transposed to other profiles. For example, for 
a Moore profile, r™ oore = 4 1 / 3 r_2, and the corresponding 
scale radius for the sub-halo is obtained from rescaling the 
NFW one, i.e. rf oorc = rf w /0.63. 

The last relation links p s to M V - 1T . Rewriting the profile 
p(r) = p s X f(r), we get 



Ps 



47T J Q v " r 2 f(r)dr 



(6) 



When sub-halos are embedded in a larger host halo, the 
virial radius does not describe the physical radius anymore, 
and the integration should be performed up to the actual 
gravitational boundary of the object. For the isothermal 
case, for instance, the bound radius would be defined 
such as p(rb) = 2ph osi (r), where r is the sub- halo location. 
Nevertheless, such a change is negligible for small clumps, 
apart from the very central regions of the host halo, of 
which all details are erased by GCR propagation. We will 
therefore neglect this further, as we have checked that it 



does not affect our results. For a NFW profile, the integra- 
tion leads to (see e.g. Fig. [2|) 



1 i'w 



M vir /{4nriA( 



:■)]• 



(7) 



where A(c vir ) = ln(l + c vir ) - c vir / (1 + c vir ). 

Hence, in these models, once the c V i r — M V i r relation is 
specified, the profile of a clump is fully determined by its 
virial mass M vir . The behaviour of r^ fw (M vir ), p" fw , as well 
as c v ir and other related quantities, are illustrated in Fig.[2j 



2.2.1. 



M vir relation: B01 and ENS01 models 



We will use the two toy mo dels B01 ([Bullock et al.ll200l[ ) 
and ENS01 ([Eke et all 120011 ) . which are based on N-body 
simulations; we refer the reader to these two papers for a 
detailed description. 

These models predict that the halo concentration 
decre ases with the halo mass (see also iNavarro et al.l 
Il997f ). Note that this behaviour has been observation- 
ally co nfirmed recently at the c luster scale (jBuote et al.l 
120071: IComerford fc Nataraianl I2007D . albeit with a 



slightly higher c v 



M v 



norma lisation than pre- 



dicted (|Comerford fe Nataraianl [20071). In subsequent No - 
body simulations (IWechsler et al.ll2q02l: iKuhlen et al.|[2005l : 
IWechsler et~aT1l2006t IMaccio et al.l l2007a|) . a good agree- 
ment was found with the B01 model, adjusted to a slightly 
lower normalisation ~ 15 — 20%Q On the other hand, the 
ENS01 model is excluded in lMaccio et all (|2007al ) (because 
of a too shallow slope), but preferred from the analysis of 
analytic lens models in iFedeli et al.l (|2007l ) . The state-of- 
the-art results for halo c oncentrations com e from the recent 
Millennium Simulation (jNeto et al.ll2007l) . The analysis at 
3 = shows a clear disagreement with B01 for high halo 
masses M vir > 10 13 M Q with a better match with ENS01. 
Still, no conclusion can be d rawn to fayour o ne model or 
another at the low mass end ([Neto et al.ll2007t) . 

In any case, both to y models are likely to be no t realis- 
tic enough. For cxample. lRomano-Diaz et all (|2007h showed 
that, to some extent, the evolution of c v i r could forget the 
initial conditions depending on the degree of violence in 
its merger events. The dependence of dark halo cluster- 
i ng on the concentration parameter also affects the relation 
(ISheth fc Tormenll2004t IWechsler et aLll2006l: iGao fc White! 
120071: Uing et al.l l2007t lDiemand et al.ll2007bD . but this is 
sub-dominant compared to the observed dispersion of c v i r 
([Ragone-Figueroa fc Plionisll2007l : IMaccio et aT]l2007al and 
references therein). 

However, as our goal is to bracket the uncertainties due 
to the clumpy contribution, we will stick to the simple de- 
scriptions of B01 and ENS01, which give respectively an 
upper limit and a lower limit on the concentration for t he 
lower masses (see also Fig. 1 in IColafrancesco et al.l[200 6^Fl. 
For our purpose, it is sufficient (and convenient) to use a 
fitted polynomial form at z = 0, to encompass the two 
extreme cases (see Fig. [2] for an illustration of the B01 re- 



4 The concentration is mod ified when consid ering dark en- 
ergy with various values of to (|Dolag et al.ll2004h . the constant 
equation of state, bu t still remains consistent with B01 model 
(IWechsler et al.ll20ul ). 

5 Their figure corresponds to slightly modified B01 and 
ENS01, which are not retained here. 
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lation): 



with 



ln(cvir) = E C * X I ln (^)j 



(8) 



Cf 01 = {4.34, -0.0384, -3.91 x 10~ 4 , -2.2 x 10~ 6 , 

-5.5 x 10~ 7 } (9) 

and 

c ensoi = | 3 14) _ 018j _ 4 06 x 10 -4 j Qj o} . (10) 

2.2.2. Scale radius and local DM density in the Milky Way 

The DM smooth halo of the Galaxy also follows the pre- 
vious relation, but a more precise description of the halo 
properties is in principle possible from in-situ observations. 
However, the determination of the mass distributions in the 
Galaxy remains challenging, especially if no assumption is 

made on the DM profile . 

In lDehnen fc Binnevl (|l998f ). the authors fitted a multi- 
parameter mass model to the available kinematic data for 
the Galaxy. They found a wide variety of models surviv- 
ing the fitting process, showing that the mass distribution 
within the Milky Way is still ill-dcfincd. The case of the 
NFW profile, which performs as good as any other profile 
(see Model 2c in their Table 4), was best fitted with a scale 
radius r s = 21 .8 kpc and pp, ~ 0.27 GeV cm -3 . A more re- 
cent analysis (iBattagli a et al.|[2005l ) using a new sample of 
240 halo objects (including globular clusters, satellite galax- 
ies and stars) found that the isothermal profile was ruled 
out for a constant velocity aniso t ropy of DM. This conclu- 
sion is disputed bv lDehnen et al.l(|2006l ). Nevertheless, both 
analyses come to similar conclusions concer ning the best fit 
xofiles. For example, for the NFW profile. IBattaglia et al.l 
20061) find r s = 17.3 kpc, p G ~ 0.4 GeV cm" 3 (corre- 
sponding to A/vir = 1-5 x 10 12 M Q for R vir = 312 kpc). 
These values are in agreement with those f ound for simu- 
lated halos with similar mass ranges (see e.g. lNavarro et all 
120041: iDiemand et al.ll2007ah : in the Via Lactea run, r s ss 
25 kpc, po ~ 0.2 GeV cm -3 (corresponding to M v ; r = 
1.77 x 10 12 M^ for R vir = 389 kp c). Finally note that a 
recent analysis ([Cowsik et al.|[2007l ) based on the dynamics 
of dwarf-spheroidals derives p Q w [0.25 — 0.4] GeV cm~ 3 . 

We shall fix the parameters of the Milky Way smooth 
DM distribution. In agreement with the previous values, we 
define our reference model with r s = 20 kpc and normalise 
it to the local density pq = 0.3 GeV cm~ 3 at i? = 8.0 kpc. 
This allows the calculation of the mass within any radius r, 
which gives the virial halo radius i?^ ir when combined with 
Eq. ©: = 280 kpc (so that c vir = 14), = 1.1 x 
10 12 Mq. Varying p Q only changes the overall normalisation 
of all fluxes, whereas modifying r s would slightly change the 
spatial distribution, which docs not affect the conclusions 
of this paper. 

2.3. Number density of clumps n c \(M c \, r) 

High resolution simulations have revealed that a large num- 
ber o f self-bound substructures survived in the Galactic 
halo dTormen et al.lll998l: iKlypin et al.lll999t [Moore et all 



[19991: iSpringeT et al.l " 120011 : iZentner fc Bullock! 120031: 



De Lucia et al. 2004; Kravtsov et"aTll2004l: IWeinberg et all 



2006; Ma ccio et al]|2007al [Dicmand ct al. 2007a) 



The mass and spatial distribution of sub-halos shown 
by these simulation s can be approximated as (e.g. 
IDiemand et~alll2004bl ) 

dN d (r,M d ) _ dn cl {r,M cl ) dV M {M cl ) dV v (r) 

— Jv c i X tt^ X 



dVdM c 



dM r 



dM r 



dV ' 

(11) 

where the last two quantities are probability functions: 



A/min 



dV M (M cl ) 
dM r , 



dM c] = 1 



dTyjr) 
dV 



dV = 1. 



(12) 



(13) 



The parameter N c \ is the total number of clumps within the 
virial radius R^ ir of the Galaxy (see Sect . 12. 3. 3")) . This means 
that the mass distribution of sub-halos does not depend of 
their locations in the host halo (tidal effects modify this 
picture, but only in the very central regions of the host 
halo ) . 

2.3.1. Mass distribution and cut-off 

For the mass distribution, the following power-law depen- 
dence is observed: 



dV M (M c i) 
dM r] 



K M x 



Md 



(14) 



The factor Km is such that the previous distribution is 



normalised to 1 for M c \ G [M TO i n ,M max \ 



KM = W G X 



(ot m - 1) 



M 



l-c„ 



In the 
DM 



limit 



Mr, 







> M n 



M u 



M ; : 



(15) 



we have K 



M 



Moore et all 



("m - 

The l ogarithmic slope q m sa 2 (e.g . 
1991 lGaoet~aTI 120041 : iReed et all 120051 : 
Diemand et al. I I2006D . but the range of a m values ob- 
tained in publish ed studies spreads between 1.7 and 2.1 
(HE aw et al. l l2007|. and references therein ). For their Milky 
Way simulation. IDiemand et all (|2007al) find - 1.9. Note 
that this is in agreement w ith the value a m = 1.91 ± 0.03 
found in IShaw et alj (|2007l ). However, when using an im- 
proved identification method of sub-halos on the same sim- 
ulations, the latter authors conclude to a shallower depen- 
dence a m = 1.79 ± 0.04. 

The mass distribution covers a wide range, from the 
heavi est sub-halo mass in the Galaxy, M max ~ 1O 1O M 
(e.g. iMoore et all Il999t IDiemand etal] l2005ah . down to 
a mass M m j n , of which the value is still debated. At an 
early stage of structure formation, a cut-offs on the lower 
masses appears due to (i) the diffusion of the DM parti- 
cles (collisi onal damping) out of a fluctuation and ( i i) free 
streaming dHofmann et al.l 120011 : iBerezinskv et all 120031 : 
iBringmann fc Hofmannl 2007 , and references therein) . The 
first process occurs after freeze-out of the DM particles, 
when it is still in kinetic equilibrium for some time with the 
thermal bath (leptons, quarks, gauge bosons). Elastic and 
inelastic scattering on fast particles results in momentum 
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exchange such that DM particles diffuse in space, leading to 
a cut-off mass Md for the structures. After kinetic decou- 
pling, the particles move freely in the expanding Universe 
background and the temperature of this decoupling sets the 
free streaming cut-off Mf s of the mass spectrum. Both cut- 
off depend on the DM cand i date propert i es. For neu trali- 
nos. lBerezinskv et al.l (|2003t) ; lGreen etaTl (|2004l . l2005f) find 
M D « 10" 12 - 10- 10 M o and M fs w 10" 8 - 1O" 6 M . This 
lower mass is slightly increased when taking into account 
acoustic oscillations owing to th e initial coupling between 
the CDM and the r adiation field (jLoeb fc Zaldarriagal2005t 
iBertschingen 120061 ). A more careful analysis of the tem- 
perature of kinetic decoupling taking into account a more 
realistic range of variations of the particle-physics mod- 
els consistent with co smological data was recently done in 
iProfumo et all (l200l . Considering SUSY models (MSSM 
and mSUGRA) as well as models with universal extra di- 
mensions (UED), these authors found the range Mf s £ 
[10~ 12 - 10" 4 ] M . 

To follow t he hist ory of these tiny substructures, 
iDiemand etafl (|2005aft performed a high resolution N- 
body simulation. The authors were able, for the first time, 
to resolve a Milky- Way size dark halo down to the free- 
streaming stage. They report survival from the small- 
est structures (injected down to M m i n ~ 1O _6 M , size 
~ 0.01 pc) at z = 26. However, tidal destruction of the 
lightest clumps and encounters with sta rs are still possi- 
ble at late stages. In an analytical model. iBerezinskv et al.l 
(|2006| ) compared the strength of tidal stripping (i) dur- 
ing the hierarchical clustering, (ii) by stars from the stellar 
bulge, (iii) by stars from the halo and (iv) by the Galactic 
disk. They found that the last of these processes was the 
most effective, predicting that only 17% of the Earth-mass 
clumps survived the tidal destruction. Note that the ef- 
ficiency of tidal disruption depends on the mass of the 
clump but also on its environment (position in the Galaxy) 
so that, in principle, Eq. (fTTj) cannot be used. Indeed, 
tidal stripping is more efficient towards t he Galactic cen- 
tre: for example, [Berezins kv et al.l (|2006t ) predict no light 
clumps at the radial distance r < 3 kpc. However, the frac- 
tion of surviving clumps is still controversial. Several re- 
cent studies have focused on the fate of these Earth-mass 
clumps. Alt hough some of them conclude to near-co mplete 
destruction (jZhao et al.ll2007t [Angus fc Zhaol I2007D some 
others underlined their resilience (iHavashi et al.l l2003t 



iGreen fc Goodwin 2007; Goerdt et all 2007^ in the Galactic 
potential. In the latter case, it is likely that the inner den- 
sity slope of cuspy satellite halo s remains unchanged, eve n 
if the halo loses a lot of its mass (jKazantzidis et aLll2004bf) . 

In any case, as we have already emphasised, the contri- 
bution of the central regions of the Galaxy is suppressed 
by the diffusive transport (for charged particles), there- 
fore it is expected to be unimportant. This assumption is 
reinforced by the fact that, compared to the smooth dis- 
tribution that is cuspy, the clump distribution might be 
cored (see Sect . 12 .3721 below) . We checked that taking or not 
taking into account a significant destruction o f low mass 
clumps — as modelled and described, in e.g. |b| 120061 — left 
the results unchanged. Thus, for our purpose, Eq. (|16|) is 
a good enough description of the clump distribution. The 
mass distribution is then fully characterised by its slope a m 
and its minimal mass cut-off M m j n . 



2.3.2. Spatial distribution of clumps 

In most N-body experiments, the spatial distribution of 
clumps is found to be anti-biased with respect to the DM 
density, at least down t o the smallest clumps resolved 



10 6 Mm) at the m o ment dGhigna et al.ll2000t 
2004 iGill etal] 12004 iGao et al.l 12004 



De Lucia et al 



Diemand et al 



2004bl . and references therein) . It is parametrised as (spher- 
ical symmetry is assumed) 



dV v (r) 
47rr 2 dr 



= K v x 



1+1 — 



(16) 



where the core radius, is a fraction of the virial halo 
radius R^ h . . The constant K y is chosen here to ensure nor- 
malisation to unity when integrating over B^ n : 



H 



nh 



tan 



nh 



IDiemand et all (|2004bh found r H ~ 0.14 R^ ir for galactic- 
like sub-halos. This bias could be due to the fact that, on 
average, tidal mass loss experienced by sub-halos is larger 
in the inner regions than near and beyond the virial radius. 
This resu lt seems to be largely una ffected b y the baryon dis- 
sipat ion l|Nagai fc Kravtsovl [20051 . but see IWeinberg et al.l 
120061 for a slightly different conclusion from a SPH simula- 
tion) . 

However, some recent studies ar gue that this cored dis- 
tribution could be a selection bias dKuhlen et alj 12007!) or 
a limitation of co llisionless simulations (iMaccio et al l 20061 : 
IShaw et al.ll2007h . For example. iKuhlen et all (|2007h find 
in their Via Lactea run a spatial distribution that matches 
the prolat e shape of the host h alo. The same trend is ob- 
served in iMaccio et al.1 ([2006), where the dissipation of 
the baryons greatly enh ances the survival of the sub-halos. 
These authors (see also iNagai fc Kravtsovll2~005h find that 
the clumps profile is well fitted by a NFW, even if the latter 
is still less concentrated (c v i r ~ 6.5) than their simulated 
overall mass distribution (c vir = 9.6). Indeed, the smallest 
clumps are likely to follow the smooth DM spatial distri- 
bution, and such an assumption has very often been used 
in analytical studies of DM dumpiness effe cts on gamma- 
ray production (e.g. iBerezinskv et al.ll2003l ). For the sake 
of completeness, such a configuration will also be used later 
for the calculations, and to be conservative, the space dis- 
tribution of clumps will be taken to be exactly that of the 
smooth component (same global concentration relation). 



2.3.3. Clump number normalisation N c \ 

The parameter N c \ is often determined by adopting the 
nu mber of sub-halos w ithin a mass range. For exam- 
ple, M oore et all (|l999l ) found 500 sub-halos with bound 
ma sses > 1O 8 M . The recent Via Lactea simulation 
of IDiemand et~aTT (|2006h gives N d (> M ref ) = 6.4 x 
10~ 3 (M re f/1.8 x IO^Mq) -1 , which corresponds to N d (> 
1O 8 M ) « 115. 

In a more general con text of various masses of host 
halos. several simulations (Ivan den Bosch et al.l |2005, and 
references therein) are compatible with the value N c \(> 
M ref ) = 0.017 x (A/ rof /A/ host )-°- 91 . Taking a mass M host = 
1.1 x 1O 12 M for the Galaxy leads to iV c i(> 10 8 A/ C i) « 81. 
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For definitiveness, we choose to set the normalisation 
N cl such as N c i(M c i > W 8 M Q = M Te{ ) = 100 = N iei . 
Taking an upper bound of M max = 10 10 M Q , we get for 
a m ^ 1: 



N ie 



1) 



Km Ml~r 



(iVrefM^" 1 ) X M D 



max 
l-a n 



(17) 

(if M min < Mnax) (18) 



where Km is the normalisation given in Eq. (|15|) . For in- 
stance, taking {M min , M rof , M max } = {10" 6 , 10 8 , 10 10 } M Q 
and a m = 1.9, we find N c \ ~ 4 x 10 14 clumps , consistent 
with values obtained bv lDiemand et all (|2005aj ). 

3. DM modelling choices and salient features 

Having discussed in detail the values, uncertainties and rel- 
evance of various parameters entering the DM distributions 
(both smooth and clumpy), we now summarise the refer- 
ence configurations used as inputs of this paper (Sect. COT) . 
Two main consequences are observed: the index of the mass 
distribution strongly affects the mass fraction of DM in 
clumps (Sect. [22]), whereas the c vu - — M v ii relation impacts 
on the luminosity profile (Sect. 13.31) . 

3.1. Reference configurations 

The distance of the Sun to the Galactic centre is fixed 
to i?© = 8.0 kpc. Whatever the clump configuration, the 
virial radius of the dark halo in the Galaxy is set to 
-Rvir h = 280 kpc, and the local DM density (smooth and 
clump altogether) to p Q = 0.3 GcV cm" 3 (Sect. |2~2~2|) . 



Clump description 


Values 


dVv{r)/dV 


Cored* or NFW 


Inner profile 


NFW* or Moore 




[1.8 - 1.9 1 - 2.0] 


Af min 


[10~ 6 1 - 1 - 10 6 ] M Q 




B01* or ENS01 



Reference configuration. 

Table 2. Description of the various configurations used in 
the paper for the sub-halo parameters. 



solely on c v i r and M v i r through Eqs. and j6]). The 
concentration c v i r depends on the virial mass M v j r , as 
provided by the B01 and ENS01 models (see Sect. EHTTl 
Eqs. Sand EE]). 
2. The clump numerical density n c i(M c \,r) is given by 
Eq. (pi]) with r H = 0.1 x R vh h = 28 kpc. N ci is 
set from the condition N ci (> 10 8 M Q ) = 100, with a 
clump mass upper boundary of 1O 1O M0. The logarith- 
mic slope of the mass distribution is a m G [1.8 — 2.0], 
and we will survey minimal clump masses starting from 
-Mmin = 10 _6 M Q . These last two parameters completely 
set the mass fraction /m of the virial mass in clumps 
(see Sect. [321 below), defined as M*j ot = f M M^ (note 
that /m does not necessarily coincide with /, see below). 

A synthetic view of the relevant parameters retained 
in this study are proposed in Tables Q] and [51 Note that 
all varying parameters come from the clumpy distribu- 
tion (spatial distribution, A/Min, ct m , inner profile and the 
cvir — -Mvir relation). The configurations, for which we will 
calculate the boost factors, are listed in Tabled 



3.1.1. The smooth component 

It is chosen as a NFW (see discussion in Sect. 12 . 1 .3[) with 
an inner radius r s = 20 kpc (c V i r = 14). In the ab- 
sence of any clump — we denote Pg m (r) the corresponding 
smooth distribution — and with the above values for p Q , 
i? vir h and r s , we recover = M s ° m = 1.1 x 10 12 M Q . 

The fraction / is usually defined as the fraction of DM 
taken from the smooth profile and redistributed into the 
clumps. The smooth contribution in this configuration is 
then Psm {r) = (1 - f)p° am (r), such that M am = (1 - f)M$ h . 

3.1.2. The sub-halo component 

If the spatial density of clumps is oc p sm (r), the redistri- 
bution of the fraction / of the DM into clumps is straight- 
forwardly written as p c \ = fPs m {r). Note that in this case, 
we have a local density of clumps fp@ and M*j ot = fM^ ir 
(such that Mlf + M Bm = M^J. We elaborate on the im- 
portant case when the two distributions arc different in the 
next Sect. 13.1.31 Otherwise, the clumps parameters are as 
follows: 

1. The inner profile of the clumps p c \{r) is taken as a 
NFW or a Moore. The saturation density is taken from 
Eq. ©, with p sat - 1O 19 M kpc" 3 for typical WIMP 
parameter^. The scale parameters r s and p s depend 

6 See discussion in Sect. 13.31 for the consequences of varying 
p S at and considering different inner profiles (e.g. a Moore inner 
profile). 



3.1.3. Defining the local fraction of DM / Q in clumps 

If the smooth (e.g. NFW) and the sub-halos (e.g. cored) 
spatial distributions arc different, the mass fraction of DM 
in sub-halos within r < R^ lt is not constant, but depends 
on the galactocentric radius r. The point is that in order 
to compute boost factors, one would naively want to sub- 
tract any fraction of DM added in the form of clumps to 
the smooth component, and compare this new setup to the 
case in which DM is only smooth. A clear definition of that 
fraction is crucial before going further. Indeed, we show 
hereafter that if not treated carefully, there is a source of 
ambiguity in the interpretation of the resulting boost fac- 
tor. 

Let us first introduce the total mass carried by the 
clumps within Ry- 1T , defined as 

M c \ ot = JV d / M '" aX dM cX M cl dV ^ cl) = iV a (M cl ). (19) 

The quantity (M c \) is the mean clump mass associated with 
the mass range [M m i n — M max ], the mass probability dis- 
tribution (Eq. fl"4"l) . and N c \ is the total number of clumps 
fEq. fTT]) . Without loss of generality, the total density profile 
of DM may be expressed as 

Ptot(r) = (1 - f) Psm (r) + M d ot (20) 

where / is a DM fraction subtracted to the smooth 
component, which wc discuss later on. The quantity 
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Cosmology 


DM (Milky Way) 


Clumps 


Smooth DM halo 


n m = 0.24 


i? vir = 280 kpc 


Global 1 ": cored, a m =1.9,M min =10- 6 M Q 


NFW (1,3,1) 


p c = 148 M© kpc" 3 


Mvir = 1.1 x 1O 12 M 


Inner profile: NFW, B01* 


r s = 20 kpc (cvir = 14) 


A vir (0) = 340 


p (i? Q ) s = 0.3 GeV cm" 3 


AT cl (> 1O 8 M ) = 100 


Psat* ~ 10 19 Af Q kpc" 3 



Number density of clumps as denned in Eq. JlH , 
c v i r — M v ir relation, see Eq. J8j. 
^ The distance of the sun to the galactic centre is set to Rq — 8.0 kpc. 
Full expression for the saturation density is given by Eq. l[2j. It also applies to the clumps inner profile. 



Table 1. Useful parameters and reference configuration of the DM modelling (see text for further details). 



M*° t dPv (r) / dV is merely the averaged mass density profile 
of the whole sub- halo population, obtained from integrating 
Eq. pip over the whole mass range of clumps. 

Two observational constraints can help to define what 
kind of fraction / is needed for consistency: the total mass 
of the Galaxy M^ h ., and the local density pq. If one wants 
to ensure that the total mass is left unchanged when adding 
clumps, then / is the mass fraction /m given by: 



/ = /. 



M 



(21) 



Otherwise, if one prefers to have a constant local matter 
density, then / is defined as a local density fraction / as 
follows: 

/ = /eS ^ x fM. (22) 



dV 



First, if the spatial distribution of clumps tracks the 
smooth profile, then we have by definition dVy(r)j dV = 
Psm(r)/M^ ir , and /m = /©■ Using either one of these frac- 
tion concepts is therefore equivalent: in other words, the 
halo mass and the mass density (at any r) are conserved 
when sub-halos are added. 

Now, if the two distributions spatially differ, it is no 
longer possible to fulfil both constraints. We have no choice 
but to abandon either the halo mass to be constant, or the 
local density to be constant. Let us see what happens when 
one of the two above conditions, Eq. (f2"Tj) or Eq. (f2"2"| . is 
plugged in Eq. (|20|. 



First option — ensuring M^ ir is constant: this choice sets 
that the fraction to be used is the mass fraction defined 
in Eq. (|2"Tj) . An immediate consequence is that the total 
local DM density is now slightly modified. Plugging this 
condition (Eq. \TI§ in Eq. fl20|) for r = R Q leads to 



Ptot(-Ro) 



A I 



1 



M v h ir dV v (R Q ) 



dV 



xp Q . (23) 



Taking our reference configuration, i.e. the NFW profile for 
the smooth component (Sect. I3~l~l"j) and a cored profile for 
the clump component (Sect. [3"T 2[) . gives 

Ptot(i?o) = [1 - fu (1 - 5.5 x 1CT 2 )] x p Q ~ (1 - f M ) PQ . 

(24) 

Such a modification of the local DM density translates the 
fact that the clumpy contribution is locally almost neg- 
ligible, in that case. This will occur every time the clump 
distribution is flatter than the parent one. As the local den- 
sity may certainly vary within a factor of two (Sect. |2.2.2|l . 
even putting up to 50% of the DM mass in clumps is accept- 
able. Such a hypothesis is very often made in the literature 



dealing with 7-rays, but the previous consequence is almost 
never mentioned. 

However, this choice is not judicious in our study. Doing 
so would even bring additional confusion to the issue of 
boost factors. Indeed, unlike 7-rays, we remind that for 
primary cosmic antimatter, the flux is very sensitive to the 
local density (see App. IA|) . Assuming for a while that the 
smooth component locally dominates the clumpy one (it 
will actually be shown later to be the case, see e.g. Fig. [5]), 
the calculated mean boost factor would be B e g ~ (1 — 
/m) 2 (see Eqs. Irl71and [rj5|) . which would result in a number 
significantly less than unity. This would consequently lead 
to a damping factor instead of a an enhancement, which 
would bring about misleading interpretations. 

Second option — ensuring p(Rsol) to be constant: to avoid 
the above situation, we may use the concept of local density 
faction instead, defined by Eq. (|22|). It comes to demand- 
ing ptot(-R©) = Pq- The boost factor now asymptotically 
goes to (1 — /©) 2 ~ 1 if the clump contribution is negligi- 
ble. However, this normalisation has to face again an un- 
avoidable issue: the total halo mass will be modified by the 
adjunction of clumps to the smooth component. Plugging 
back /q in Eq. (120p , integrating over the virial volume and 
using again fu (see Eq. [21]) , we get 

M tot = (1 _ Jq + Jm) x M h r ^ (1 + j m) x M h r _ (25) 

With the different clump configurations used throughout 
the paper, the total halo mass within the virial radius can 
be increased up to fu ^ 50% level (see curves a m < 2.0 in 
Fig. [1]). Such values remain within current estimates of the 
total mass of the dark halo, as recalled in Sect. I2.2."2l 



Closing the case: we conclude by reminding that neither 
of these choices is better than the other. Both are, some- 
how, equally artificial. Indeed, there is only, if so, one true 
distribution of smooth and clump DM in the Galaxy. The 
ambiguity appears because we wish to compare the calcu- 
lated fluxes to a hypothetical configuration with no clumps. 
In the context of antimatter fluxes, as explained, the sec- 
ond option (ensuring the same local DM density whatever 
the configurations) makes more sense, as it leads to boost 
factor values asymptotically reaching 1. This second option 
is retained throughout the papefl 

Finally, before closing the DM section, let us discuss 
how the various configurations gathered in Table [5] impact 



7 Let us say it again: should the first option have been re- 
tained, we would have ended up with boost factors smaller than 
one! 



10 



Lavalle et al.: Antimatter GCRs from DM annihilation: Abandoning hope in dumpiness enhancement? 



=fii'ini| i 



</> 
re 

E 

3 

E 

3 

o 



10" 1 



................ 



iiiiii| i iimiij i iiimij i iiiii^ 1 1 imiij i iiimij 1 1 mi i niiiij i ™ 

dP M / dM °= M" a - 
-°-a = 2.10 
"a = 2.05 
a = 2.00 
— a = 1.95 
♦a = 1.90 
S \ a = 1 .80 

\ \ ■■>■■ a = 1.70 

\ h 

*v V s . 



- +-+ +-+ 4— I- 4-+ 4 





i. 4-+ 4-+ 4-+ 4-+ 4-+ 4-+-S-'? 



I I 



J I I i.'J^iniJ 



10" 8 10" 6 10" 4 10" 2 1 10 2 10 4 10 6 10 s 10 10 

[solar unit] 

Fig. 1. Mass fraction /m as a function of M m i n , for different 
logarithmic slopes a m of the mass distribution (from 2.1 down 
to 1.7 — top to bottom curves). 



on some generic properties for the clumps (mass fraction 
/a/ and luminosities). 

3.2. Mass fraction fuj in sub-halos 

The minimal mass M m - m of the clumps able to form — and 
to survive tidal disruption — is a crucial parameter (see also 
next subsection). Along with the slope a m appearing in 
the mass distribution Eq. it sets the fraction of DM 

in clumps, f M = M^JM^. 

The evolution of with a m and Af m j n is shown in 
Fig. [TJ The behaviours are in agreement with the figures 
discussed in several simulatio ns. For example, t aking a re- 
solved mass Mmin > 10~ 6 Af g . IShaw et al.ll2007l (and refer- 
ences th erein) found / ~ 5 — 1 0% for a m = 1.8. As ar- 
gued in iDiemand et al.1 (|2007a|) . where a larger value of 
a m ~ 2.0 is preferred, 10% might be only a lower limit 
and this fraction could reach /m ~ 50%. In the extreme 
case of a slope a m = 2.1, all the DM could be distributed 
in clumps, even forbidding the existence of clumps of mass 
smaller than < 10 _2 M©. This latter configuration is proba- 
bly not realistic, so we choose to limit the study to the range 
a m G [1.8 — 2.0]. Consequently, as observed from Fig.[T] the 
fraction of mass in clumps M^ t , corresponding also to the 



additional mass added to M^ r (see Eq. 
range - 10 - 40 %. 



will lie in the 



3.3. Luminosity: a closer look on the astrophysical term 

Before plugging the propagation, it is interesting to have a 
look at the luminosity of the source terms in the various 
configurations. The total luminosity of DM sources can be 
separated into a particle physics term times an astrophysi- 
cal term: 



is 



>(E,r) = (dN/dE) x £ astro (r). 



The particle physics term is factored out by normalising to 
the local luminosity Lq = L sourcc (Ro) oc p© 2 . The relative 
astrophysical luminosity is then defined as 

C(r) = L somcc (r) / Lq — £ a stro(f)/£astro(-Ro). 

For short, below, we will continue to call this quantity the 
luminosity, and use M instead of M c \ for the mass clump. 



Smooth component: it is straightforwardly written as 



£ sm (r) 



Pa 2 



(26) 



Substructures: it is convenient to separate the total clump 
contribution as the sum of the contributions of each decade 
of mass (throughout the paper, the logarithm bins of mass 
are denoted Mi = [10* — 1O 4+1 ]M with i an integer). 
Taking into account the space and mass distributions of 
the clumps, the total and per logarithmic mass bins lumi- 
nosities are defined as: 



ptot 



^cl: 



dVy 

dV 



(r) x 



Mi 



dC c \ 
dhxM 



(M) d{ln M) 



with 



dC c \ 
dhxM 



(M) 



N cl xMx ^T-(M) x £(M) 



dM 



e(M) 





>ci(r,M)" 


L 


. pq . 



(27) 



(28) 



(29) 



(30) 



The above Eq. |29|) defines the luminosity mass profile while 
Eq. (|30[) defines an intrinsic effective annihilation volume 
for a clump of mass M. Note that introducing the mean 
value (C)m of £ over the i th mass bin, allows the luminosity 
per logarithmic mass Eq. (|28[) to be recast under a form 
where the dimensions appear more explicitly: 



£dW^clX(f) M x5W 



M 



10 i+1 M, 



lO^Af© 



dV 

dV*r 

dM£(M) 



(M) 



dM 



(31) 



(32) 



Last, it is useful to introduce a dimcnsionless intrin- 
sic boost factor B C \(M) for a clump (not to be mistaken 
with the global boost factor defined in Eq. [57]) . The latter 
compares the annihilation rate of the clump, to the rate 
that would be obtained for a clump that had the same vol- 
ume, but with a constant DM density p Q . The local intrinsic 
boost factor can be expressed as 



B c1 (M)e0)x| 



(33) 



We emphasise that such a quantity is meaningful since i) 
antimatter fluxes mostly depend on the local DM density 
and ii) unlike 7-rays, we do not look in one specific direc- 
tion, but rather integrate on the whole clump signal. 

Before concluding on resulting luminosities C sm (r) and 
C cl , let us further detail the various terms appearing in the 
clump luminosity term. 
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Table 3. Sub-halo parameters (reference configuration, i.e. in- 
ner NFW and BOl for c v i r — A/ V i r ) for different masses: virial 
radius R V1I (Eq.[3]), scale radius r° fw (Eq.[5]), scale density p° fw 
(Eq. 0, concentration parameter c v i r S01 (Eq. [5]), effective vol- 



ume £ nfw > B01 (Eq. EU 
See text for details. 



intrinsic local boost _B" 



(Eq.|33j 



3.3.1. Annihilation volume £(M) 

This quantity is a function of the mass clump M, and it 
depends on the inner profile (NFW or Moore), the c vlr — 
M v i r relation (BOl or ENS01) and the saturation density 

Psat- 



Reference configuration (NFW inner profile): for an NFW 

profile, the annihilation volume has a simple analytical ex- 
pression 



£ ntw (M) 



4tt 



(r nfw )3 



nfw 



p,;. : 



(34) 

where the scale radius r" fw and the density at scale radius 
p" fw depend on the clump mass M, as given in Eqs. ([5]) and 
0. The function ?y nfw (r) is denned as 



7? nfw (r) = 



1 



~nfw 



(M) 



where the saturation density for a NFW is given by 



„nfw 



(M) 



nfw 



(M) 



- nfw (M) 

Psat 



(35) 



(36) 



It is easily checked that £ (M) is largely insensitive to 
the exact value of p sat (~ 10 19 M Q kpc"" 3 ). 

The local intrinsic boost factor B^ W (M) can also be 
analytically expressed in terms of the virial parameters: 



B^(M) = 



M 



12np Q R 



' (3 + C vir (3 + Cvir)) 



(1 + C vir )[c vir - (1 + Cvir) ln(l + c v 



(37) 



As c v ir is only very slightly mass dependent and i? v ir cx 
M 1 / 3 , the intrinsic boost factor, B c \, is almost constant 
over a wide range of sub-halo masses. More precisely, for 
the NFW case in the BOl model, we find it to scale with 
the concentration parameter like 



1.29 x 10" 



(Cvir™ 1 ) 572 . (38) 



For illustration purpose, typical values for all above 
quantities (for various masses of clumps) are gathered in 
Tablc[3J After applying a suitable renormalisation, the same 
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Fig. 2. Same as in TableO Each parameter is normalised with 
respect to its maximum value (in the mass range displayed). 



quantities — some of which having power law dependencies 
with M — are also displayed in Fig. [2] We check that £ nfw 
roughly scales like M, a very common feature already em- 
phasised in the literature that focuses on gamma-rays (see 
the consequences in the next subsection). For the intrinsic 
local boost factor, we read off the last column in Table [3] 
that only clumps below M < 10 6 M Q may significantly over- 
come the local annihilation signal. The first conclusion that 
can be drawn is that massive clumps, unless close to the so- 
lar neighbourhood (which is very unlikely), will not be able 
to boost the antimatter signals. More importantly, the very 
same parameter allows an upper limit to the total boost 
expected to be set, b^ w < B01 < 20. 

Moore vs NFW inner profile — ENS01 vs BOl: these conclu- 
sions are very easily extended to other configurations. We 
actually find a very simple rescaling factor linking the an- 
nihilation volume £ moorc (A/) to the above £ nfw (M), when 

ref 
Psat 



10 1J M Q kpc" 



r oorc (M, Psat )~8xr tw (M) 



(39) 



This docs not depend on the c v i r — M V i r relation, and to a 
good extent, neither on the WIMP mass and annihilation 
cross section, nor on the clump mass, as illustrated in Fig. [3] 
Note that, as expected, the ENS01 configuration is lower 
than BOl, implying, for the intrinsic boost factor, the hier- 



, D nfw,ENS01 

arcny B , 



< 



7->nfw.B01 7->iiiooic.n/iN oui , t-jij 
B c\ - B c\ < U c\ 

(roughly corresponding to 0.1:1:10). Hence the former con- 
figuration will provide the minimal boost, the reference con- 
figuration the median boost and the last configuration the 
maximal boost. 



,mooro.ENS01 



imoorc.BOl 



3.3.2. dCd/dlnM 

The differential luminosity Eq. (|29|) is evaluated taking into 
account the mass distribution of clumps Eq. (fT^|) . The re- 
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Fig. 3. Ratio of the annihilation volume for three benchmark 
sub-halo models to the reference one, defined as = £ x /£ rc = 
B x /B™ . The reference configuration is an inner NFW profile 
with B01 concentration. The three benchmark configurations are 
X = (Moore, B01) (upper curve), X = (Moore, ENS01) (middle 
curve) and X = (NFW, ENS01) (lower curve). 
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Fig. 4. Differential luminosity (true units in kpc ) of the popu- 
lation of clumps dC c \/d\riM as defined in Eq. (|29|l . The curves 
correspond to the NFW-B01 reference configuration, for vari- 
ous values of a m . The curves for any other configuration may 
be obtained by multiplying these curves to the ratios shown in 
Fig. E] 



suit is shown in Fig. Q] for different values of the slope 
a m . We recover t he various trends seen in the literature 
(see, e.g. Fig. 8 of lDiemand et al.ll2007al iBerezinskv et alj 
l2003t ). In particular, the value a m = 1.9, favoured in sim- 
ulations, shows a roughly constant luminosity per decade. 
For smaller (respectively greater) value of a m , the lumi- 
nosity will be dominated by the heaviest (lightest) clumps. 
In that case, based upon the understanding gained from 
the previous discussions, the boost is expected to be small 
(close to unity, from the intrinsic boost factor). At the same 
time, the variance of the clumpy signal is expected to be 
large (light clumps add no contribution and heavy clumps 
are scarce). For larger a m , the mass M min of the lightest 
clump is crucial, because the latter drives the total lumi- 
nosity. These large a m configurations are expected to give 
the largest boost factors. The last step is to put together 
the smooth and clumpy luminosities. 

3.3.3. Luminosity profiles C sm (r) and £ l cl (r) in the Galaxy 

The last hint at small boost factors for the case of antimat- 
ter DM is given when comparing the smooth and clump 
luminosities. This is first shown for the reference config- 
uration i n Fig. [5J (top left pan el). As already emphasised 
(see also iDiemand et ahll2007a !h for the reference configu- 
ration a m — 1.9, the contributions to the average annihi- 
lation fluxes of any decade mass range will be almost the 
same on the who le range of mass clumps (see also Fig. 1 of 
I Yuan fc Bill2007f ). Around r = Rq, the luminosity is com- 
pletely dominated by the smooth contribution (~ 100 times 
more than the total clump luminosity), so that for this con- 
figuration, we may predict beforehand (i) no boost factor 
and (ii) a small variance on this boost factor. 



The logarithmic slope of the mass distribution cv m re- 
verses hierarchy in the mass contribution: for a m = 1.8 (see 
Fig. El top right), the more massive the population of clump 
the more luminous it is, whereas for a m — 2 (sec Fig. [51 bot- 
tom left), the less massive, the more luminous. The trade 
off is reached close to a m = 1.9. However, the total clump 
luminosity never reaches the level of the smooth one! In 
the best case (a m = 2), it is 10 times smaller. Nevertheless, 
taking larger a m values naturally leads to larger boost, in a 
more general context. 

Now, if we now assume that, instead of having a cored 
profile, the spatial distribution follows the parent one 
(NFW), we see in Fig. [5] (bottom right) that the situa- 
tion is more favourable for the boost factors. Keep in mind 
that this is an upper limit since the the averaged radial 
mass density profile of clumps is believed to be a flatter 
distribution than the smooth one (see Sect. I2.3.2|) . 

Anticipating the results of using ENS01 instead of B01, 
or using a Moore inner profile instead of a NFW, especially 
from Fig. [3j we already know that ENS01 will only further 
decrease the total clump luminosity (roughly by a factor of 
10 compared to B01), whereas Moore will increase the total 
clump luminosity (roughly by a factor of 10 compared to 
NFW). 

To summarise, from the general study of the luminosity, 
we might already conclude that no configuration of DM will 
lead to huge boost factors. Pushing all the parameters for 
the maximum effect, i.e. Moore inner profile, NFW spatial 
distribution, B01 and a m = 2.0, would possibly lead to a 
boost factor of a few, but certainly not a hundred. We can 
expect all other configurations to end up with a boost factor 
close to unity. The rest of the paper is devoted to the full 
calculations to confirm these expectations. 
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Fig. 5. Relative luminosity profiles as functions of the galactocentric radius r, in units of local luminosity Lq. Luminosities are 
plotted for the smooth DM contribution and for clumps in logarithmic mass bins of 3-decade width. Top left: reference configuration 
(ref). Top right: as ref but a m = 1.8. Bottom left: as ref but a m — 2.0. Bottom right: as ref but for a spatial distribution of clumps 
oc p B m{r). 



4. Propagation model 

In the Galaxy, a charged particle travelling from its source 
to the solar neighbourhood is affected by several processes. 
The scattering off random magnetic fields leads to spatial 
and energy diffusion (reacceleration) and particles may also 
be spatially convected away by the galactic wind (which 
induces adiabatic losses). 

In this paper, the f ramework used is the following (e.g. 
Berezins kii et al.lll990t ): for the transport processes we take 
a spatial independent diffusion coefficient K(E) = [3KqR. & 
(where 1Z — pc/Ze is the rigidity) and a constant wind 
V c directed outwards along z. Cosmic rays are confined 



within a diffusive halo L, such as the differential density, 
dN/dE = N, is bound by N(z = L,r) = 0. The free pa- 
rameters of the model are the halo size L of the Galaxy, the 
normalisation of the diffusion coefficient Kq and its slope 5, 
and the constant galactic wind V c (see Sect. l4.3[ ). Other pro- 
cesses (such as continuous and catastrophic gain/losses) are 
more species-dependent. Hence, although all charged parti- 
cles are propagated in the same framework, due to this de- 
pendence, the phenomenology of propagation is completely 
different for p and positrons. 

The reader is referred to iMaurin et al.l (|200T| ) for 
a more detailed presentation and motivation of the 
framework. Note that this model has been repeat- 
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edly and consistently used in seve ral studies to con- 
strain the propagation p arameters (|Maurin et al.l 120011 
l2002t iDonato et al l 12002ft and examine the cons equences 
(|Taillet fc Maurid 120031: iMaurin fc Tailletl 120031 ) for the 



standard p flux (jDonato et al 
d fluxes (|Maurin et al 



200lh 



2004, 2006a: 



the exotic p and 
Donato et all l2004t 



iBarrau et a ll 12002). 120051; iBringmann fc Salatil I2007D but 
also for positrons ( Lavalle et al.ll2007l ; iBrun et alj 20071 ). 



4.1. Propagator and flux for anti-protons 



It was shown in iMaurin et al.l (|2006af ) that neglecting all 
energy redistribution terms (energy losses, reacceleration 
and tertiary source term) provides a correct description at 
sufficiently high energy, while remaining good enough down 
to ~ GcV IS energies (better than 50% depending on the 
propagation parameters considered). This approximation is 
retained here. The only catastrophic losses for anti-protons 
are spallations — the particle does not survive the interac- 
tion. 

Denoting r to t = Sism n iSM-v.crf SM the destruction rate 
of p in the thin gaseous disk (tiism =H, He), the transport 
equation for a point source, defining the propagator, reads 
(jMaurin et al.ll2006al) : 



-KA + V c -^- + 2hT tot S(z) 



(40) 



For simplicity, we only consider a flux detected at solar 
position r© = (&© = i£©,y© = 0, z© = 0). For a point 
source S at rg, the corresponding flux only depends, in 
cylindrical coordinates, on the relative distances r = |r© — 
rs| and z — zg. The propagator <?©(?", z) is given by 



5i(r,z) 



exp 



-k v 



■ x (41) 

OC 

^ c~ 1 K (r\fk% + fef) sin[fc„i] sin[fc n (i - z)\ 

n=0 

where Kq is the modified Besscl function of the second kind. 
The quantity k n is the solution of 



2/lj^ cos kjiL — k d sm ky^L , 



and 



Wc also have 



c n = 1 - 



sin k n L cos k n L 
kn.L 



(42) 



k v = V C /(2K) and k d = 2h T tot /K + 2k v . 

For a source term g©(r, z,ff)Q{E) (origin is taken to 
coincide with solar location) the equilibrium spectrum at 
solar position is finally given by 



- x2ldzlrdr Q*{r, z) I d8 q Q (r, z, 9) 



in 



(43) 



4.2. Positrons 



Contrarily to nuclear species, there are no catastrophic 
losses for positrons. A more crucial point is that prop- 
agation of positrons is dominated by energy losses (e.g. 



iMoskalenko fc Strondll998| ). In that case, a monochromatic 
line at the source leads to a spectrum once propagated. This 
is at variance to p whose propagator for exotic sources is at 
constant energy. 

The diffusion equation that characterises the evolution 
of the positron number density N per unit energy, with a 
source term q(r)Q(E), reads 

-^(0AN + ^{§N)^ q (r )Q (E). (44) 

The first term is simply the diffusion coefficient written as 
K(TZ) fts K (E / E ) s . For simplicity, we have also neglected 
the effect of the Galactic wind (see next section for a dis- 
cussion). 

We proceed as in lLavalle et al.l (|2007l and see references 
therein). The synchrotron and inverse Compton losses can 
be written as dE/dt(E) = -E 2 /{E t e ), with E = 1 GcV 
and te ~ 10 16 s. Defining a pseudo-time 

; _ (E/Ep)*- 1 
t = TE 1-8 
and applying the following rescaling, 

N=(E/E ) 2 N and Q{E) = (E/E ) 2 ~ S Q(E), 
the diffusion equation can be rewritten as 

^- N - K AN = q(v)Q(E). (45) 
ot 

Thus, instead of finding the solution of Eq. (|44[) . we are left 
to solve the well-known time-dependant diffusion equation 
Eq. (35) ■ 

It proves convenient to separate diffusion along the ra- 
dial and vertical direction. Considering a source located at 
(x,y,z,t) detected at (i?©, 0, 0, to), the corresponding flux 
depends only on the radial relative distance r = \r$ — ro|, 
the distance of the source from the plane z = zs and the 
relative pseudo-time f = Ie — i_E s • The Green function 
£/©(r, z, f) of Eq. (|43|) is then given by: 

„2 



Go(r,z,T) 



'(f) 



exp 



4A' f 



x G 1d (z,t). (46) 



The effect of boundaries along z = ±L appears in 
Q 1d (z,t) only. For convergen ce properties, two distinct 
regimes are worth considering ([Lavalle et al.ll2007h : 

1. for sources close to us, it is best to use the so- called 
electrical image formula (e.g. lBaltz fc Edsiolll999l) : 



+ 00 



<^(z,f)=£(-ir 



Hr) 



exp 



(z„ - z) 
4AT f 



where z n = 2Ln + (—1)™ z\ 
2. for far away sources, a more suitable expression is 



(47) 



V n (0R(z) 



(48) 



where 

4> n (z) = sin [k n (L - \z\ 



I) I (even) 



4>' n (z) = sin [k' n (L - z)] ; k' n = n— (odd). 
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5 


K (kpc^ Myr" 1 ) 


L (kpc) 


V c (km s^ 1 ) 


max 


0.46 


0.0765 


15 


5.0 


med 


0.70 


0.0112 


4 


12.0 


min 


0.85 


0.0016 


1 


13.5 



Table 4. Propagation parameters giving the maximal, me- 
dian and minimal antiparticle DM fluxes compatible with 
B/C analysis. 



Coming back to the non-hat quantities, the propagator 
for a monochromatic point source is related to Eq. (jlS)) by 

S Q e+ (r, z,E *— E S ) = ^ x G e (r, z, f = t E - t Es ) . (49) 

It follows that for a spatial and spectral distribution of 
sources qQ(r)Q(E) (origin of coordinates at solar neigh- 
bourhood) , the equilibrium spectrum at solar position and 
energy E is given by 

f'L />oo 

(E) = — x 2 / dz rdr (50) 
4tt Jo J 

poo p2tt 

/ dE s [Q(E s )ge + (r,z,E<-E s )\ / d6q{r,z,e) . 
4.3. Propagation parameters 

A few important points are reminded concerning the role of 
the various transport parameters on the propagated spectra 
of the antiparti cles created in the DM halo. More detail s 
can be found in lDonato et all (|2004l ); iBarrau et all (|2005f) . 

The halo height L determines the total number of 
sources inside the diffusive region and the typical distance 
a GCR can travel before escaping from the Galaxy (see also 
App. [S]). The galactic wind wipes the particles away from 
the disk, and a similar effect occurs if V c is large enough. 
The parameters L, V c and Kq are correlated. In the subset 
of paramet ers giving the observed B/C ratio (|Maurin et al.1 
120011 12002D , low values of Kq generally correspond to low L 
and V c , so that the DM signal is expected to decrease with 
decreasing Kq. On that basis, extreme and median param- 
eters can be extracted, in the sense that these parameters 
lead to the minimal and maximum expected flux, while the 
median parameters (best fit to B/C data) provide the most 
likely flux. These parameters are recalled in Table 0] 

Having in mind the connection between the propaga- 
tion parameters and the fluxes, we can now justify dis- 
carding, for our calculations, the effects of the wind and 
reacceleration for the positrons. For example, for configu- 
rations with small 6, as the effect of the wind is always 
negligible for anti-protons, it is also the case for positrons 
(their travel time in the Galaxy is less or at most that of 
the anti-protons). For the sets of parameters with larger S, 
the effect of the wind becomes dominant below < 1 GeV. 
However, we are mainly interested in the high energy regime 
for positrons. Furthermore, if the low-energy behaviour is 
strongly dominated by convection (as is the case for anti- 
protons when (5 = 0.85), then it superseeds energy loss ef- 
fects for positrons: in that case, all the conclusion about p 
would also hold for e + . 



5. Methods 

The smooth contribution is straightforwardly calculated, 
contrarily to the clumpy contribution that is plagued by 
statistical uncertainties (in the sense that the position of 
clumps is a random variable, see Sect. 15. 2[ ). The latter issue 
is the primary concern of this section. 

Two complementary approaches are followed to calcu- 
late the Galactic variance of the clumpy contribution. The 
first one (|Lavalle et al.ll2007t) is a semi-analytical calcula- 
tion of the mean and variance from the generic statistical 
properties of the clumps (spatial and mass distributions), 
using the particle propagators that we recalled. The second 
one uses the same ingredients, but quantities under scrutiny 
(mean and variance) are obtained by accumulating realisa- 
tions of a clumpy galactic halo. Due to the lack of any clue 
about the precise location and intrinsic properties of each 
individual DM clump, working with statistical tools is well 
motivated. The numerous clumps can be treated as ran- 
dom objects, which average properties are taken here from 
N-body simulations. 

Note that both methods rest on the assumption that 
clumps arc considered as point-like sources. This is correct 
while the distance of a clump to the Earth is greater than 
its spatial extension, and if the GCR propagation proper- 
ties do not change within the spatial extension of a clump 
(see Table [3]). As the flux, on average, is not dominated by 
nearby substructures, and since for those far away clumps 
the spatial dependence of the propagator is smooth enough 
(diffusive process), the point-like source assumption holds. 
Would a nearby clump dominate the positron or anti-proton 
flux — which is very unlikely according to our calculation — 
a single source computation would be enough to deal with 
the dumpiness issue. Nevertheless, such a case, while easier 
to calculate, would make the dumpiness itself an absolutely 
unprcdictive scenario for the indirect search for DM using 
antimatter GCRs, and is beyond the scope of this paper. 

5.1. Generalities 

Before exposing the methods, it is convenient to define a 
pseudo-Green function, denoted Q, by absorbing the energy 
dependence of the GCR propagators. To this aim, we de- 
fine the quantity dN/dEs(Es) to be the antimatter species 
spectrum at the source, which is defined here as the num- 
ber of antimatter particles injected per annihilation and per 
energy unit. 

For anti-protons, the pseudo-Green function reads: 
d,N 

g P (E)^—(E)xgg(r,z,E). (51) 

The propagation term and the source term can be factorised 
(no energy mixing during propagation, a p emitted at E$ 
is detected at the same energy E$). It means that the re- 
sults for the relative uncertainties on the fluxes and for 
boost factors arc independent of the particle physics model. 
Unfortunately, this is not the case for positrons: 

pOO f J AT ^ 

G e+ (E)^ J e dE s } [ — (E s )xg^(r,z,E^E s )j . 

(52) 

The integral characterises energy losses and the source spec- 
trum cannot be factored out of the integrand. Nevertheless, 
in the following, in order to keep the discussion at the 
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most general possible level, we will mainly focus on a 
monochromatic line of positrons at Es, i.e. dN / dEs{Es) = 
S(E — Es). The results for positrons are thus forced to be 
independent of any particle physics model. As we will sec in 
Sect.[SJ as the boost factor is close to unity for all energies, 
convolving the propagator with a realistic DM source spec- 
trum would yield a similar boost factor as obtained from 
the monochromatic line. 

In the following, we will make use of Q 1 where the energy 
dependence is implicit for any species. 

The total GCR flux cj> tot originating from DM annihi- 
lations may be separated into two contributions, for the 
smooth component and for clumps: 



"tot — <^sm 



Smooth contribution: The smooth contribution to the flux 
(j) sul is calculated using the smooth density profile p sm : 

2 

x) , (54) 



f S 
Air 



d 3 x 



Pq 



where v is the cosmic ray velocity, Q is the pseudo-Green 
function, defined above for p and e + , and S is a particle 
physics coefficient^ depending on the WIMP modeo 



S = 5- 



Pa 



t-wimp 



(55) 



The normalisation is thus chosen with respect to the local 
DM density p Q = 0.3 GeV cm" 3 and 6 = 1 (1/2) if the 
WIMP of mass m w i mp is a Majorana (respectively Dirac) 
particle. S is actually counting the number of annihilations 
occurring in an infinitely small volume, in which the DM 
density is set to p®. 



Sub-halo contribution: The Galactic halo is populated by a 
constellation of many clumps, whose positions and masses 
are actually unknown. Nevertheless, if N c \ is the number 
of clumps in a certain diffusive volume in the Galaxy, then 
their total contribution to the flux reads: 



A',i 

£< 



(56) 



The spatial dependence of the propagator is smooth enough 
(diffusive process) so that Q may be considered constant 
over the clump scale: each clump behaves as a point-like 
source. The cosmic ray flux measured at the Earth from 
the i th clump is therefore given by: 



i(x©) = T~ x S x & x ^( x © 



(57) 



where is the effective annihilation volume defined by 
Eq. O. 

5.2. Semi-analytical calculation of the flux and the boost 
factor, and associated variances, due to sub-halos 

In this secti o n, w e apply the formalism developed in 
lLavalle et al.l (|2007t ) in order to predict how boosted the 
antimatter cosmic ray fluxes should be when adding sub- 
halos. 



8 So that dN/dE x S = Q(E), as used in Sect. H 

9 This differs from the convention used in lLavalle et al.l (|2007f ) 



for which (v/in) is included in S. 



5.2.1. Whole sub-halo flux 



The propagator describes the probability for a cosmic ray 
injected at position xs with energy Es to be detected at the 
Earth (x Q ) with energy E (recalling that for anti-protons, 
Es = E, as they do not loose energy). 

As the intrinsic luminosity of a clump is entirely set once 
its mass is known, the effective volume £j can be expressed 
as ^i(M c i ; i). Thus, given Eq. (|57|) . the flux associated with 
a single clump is a stochastic variable that depends on two 
probability distributions: the space and the mass distribu- 
tions (Sect. [2~3"lh This is summarised in the following equa- 
tion: 



(53) 2l± 



dPcf, dP(x c i,M c ; 



1 



d 3 x c \dM c 



dV v (r) x dV M (M cl ) 



Airr 2 dr 



(58) 



where both distributions, given by E ns. (fl6| ) and (fl4|) . re- 
spectively, are considered uncorrelatccPi. 

The halo is populated by a constellation of many clumps 
whose total contribution to the GCR flux is given by 



0cf(x ) 



N c 

?:=i 



V 
4-7T 



x^(x 



© 



(59) 



i=i 



Though the previous expression would be the actual ex- 
pected flux for our Galaxy, we do not neither know the num- 
ber nor the precise locations and the masses of clumps in 
the halo. Nevertheless, the knowledge of their phase space 
distribution can be used to determine the mean value of 
that flux: 

N cl x ^- x Sx (£> M x (g) v , (60) 



*cf> 



47T 



where AT c i is the number of sub-halos hovering in the DM 
volume V of interest, and 



(G)v = (9) 



d 3 x £(x Q <- x) x 



(Om = (0 = 



Mr, 



M„ 



dM £(M) 



dPy(x) 
dV 

dV M (M) 
dM 



(61) 



(62) 



Equation (|60[) is the mean value of the flux, in the sta- 
tistical sense, due to all clumps in the Galaxy for a given 
model (space and mass distributions f^l. Anticipating the 
next section, we stress that any MC approach should con- 
verge to these values when taking a very large number of 
halo realisations. 



5.2.2. Variance cr*"* of the whole sub-halo flux 

The fact that we do not know how clumps are actually 
distributed, in the phase space defined by their locations 
and masses, can be expressed in terms of a variance tr*"* 
associated with their total mean flux (0*°*). For a single 
clump, the relative flux variance is given by: 



o~ 



2 



0~ 



(T 



x CT « 
2 X 7^2 



(G) 2 (0 



We remind that tidal disruption of a clump in the Galactic 
centre depends either on its mass and on its location, which in- 
duces a small correlation between the mass and the space distri- 
butions. Nevertheless, we have checked that it could be neglected 
for this purpose. 

11 The integration volume V is the DM halo volume, but in 
practical calculations, we reduce it to the diffusion volume. 
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where the individual variances affecting Q and £ are respec- 
tively 

^^(x^xjx^ - (G) 2 ; (64) 



halo 



and 



M., 



dV 



dV M 



dM i 2 {M) x 



(0 2 • (65) 



Quantities related to £ and Q will be quoted as mass- and 
space- related, respectively (see Sect. [6]). We will further 
show that mass-related effects dominate among the rela- 
tive variances, in spite of sizable space-related ones. The 
(third) crossing term of the right hand side of Eq. ([63)1 will 
consequently set the full variance of the clump flux. 

The resulting relative flux variance for the whole popu- 
lation of sub-halos is then merely: 



^tot 

4® 



1 Ccl 



(66) 



5.2.3. Boost factor B c r and its variance o~b 



Once the contribution of sub-halos to the flux is fully 
determined, the boost factor is easily computed for 
any species. As cosmic ray propagation has an ex- 
plicit energ y dependence, the b oost factor is also energy- 
dependent ([Lavalle et al.l |2007[ ), and, of course, also d e- 
pends on the cosmic ray species ([Maurin fc Tailletll2003[ ). 

The energy-dependent mean effective boost factor is 
given by the sum of the clumpy and the smooth contribu- 
tion divided by the flux that would provide the only smooth 
reference halo Pg m (r) (see Sect. 13-lTTj) : 



tot 

B eS = (l-f Q f + ^ 



(67) 



where the local density fraction / Q has been defined in 
Eq. (|2"2"|) . in order to keep the local DM density constant 
(see Sect. EH}. 

It may be useful to determine the limit for which only 
an infinitely small volume around the Earth 5(x — x Q ) is 
taken into account. Actually, this will give a rough estimate 
of the asymptotic (maximum) value of the boost factor for 
both positrons (at detected energies very close to injected 
energies) and anti-protons (at low energies), because we are 
blind to contributions from regions close to the Galactic 
centre — where the smooth DM density dominates — in 
this case. This local asymptotic value is given by: 

dV 

B e = (1 - / ) 2 + N cl x (Om x jy(R®) • (68) 

This expression neither depends on the WIMP model, nor 
on the species. We see that, as dV/dV(Ro) = 3.9 x 
10~ 7 kpc~ 3 in our reference model, only configurations with 
N c i x (£)m ?J 2.6 x 10 6 kpc 3 will yield a relevant averaged 
contribution of clumps compared to the smooth component. 
From Table [31 one can already see that the corresponding 
probability is likely to be very small. We can therefore pro- 
vide a very simple criterion to check whether a boost is 
likely to appear in any (many-object) configuration: 



n c i(R & ) x (0 



N^(Rq) x (0 > 1 
=> B cS > 1. 



(69) 



Should the sub-halos spatially track the smooth compo- 
nent, then one would get: 



S»° = (1 - f M ) 2 + hi x 



(M c l) 



(70) 



where Jm is the mass fraction of DM in clumps, and 
(M c \) is the mean mass of clumps. The previous crite- 
rion to get a relevant boost factor then becomes merely 

(0 > (Mcl)/(fM P0 ). 

The above value of the boost factor fluctuates up to a 
variance erg, which reads 
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'cl 



leading to 



a 



lot 
cl 



B ~ (1 - /o ) 2 </>! 



6 tot 



(71) 



(72) 



If the sub-halo contribution dominates over the smooth 
component, then the relative variance of the effective 
boost factor is roughly equal to that of the sub-halo flux. 
Nevertheless, as soon as sub-halos become irrelevant in the 
flux estimate, the variance of the boost factor is strongly di- 
luted by the smooth term. In this case, we obviously find a 
very small variance associated with the boost factor, even 
when the relative statistical uncertainty on the sub-halo 
flux itself is large. 

5.3. MC approach 

A complementary approach is to calculate and add ex- 
plicitly the contribution of each clump by MC drawing. 
Simulating many realisations of the DM sub-halos is an- 
other way to extract the mean flux as well as the variance 
of the clump contribution. The ensuing calculation of boost 
factors is as before, but in addition, MC provides the law of 
probability for the stochastic variable (j> c \ that describes the 
single clump flux, which is hardly inferred from the clump 
phase space distribution itself due to the needed convolu- 
tion with propagation. 

From a technical point of view, it is very inefficient to 
calculate contributions from so many sub-halos (e.g. > 10 15 
for the lightest ones) one at a time. Indeed, for the clumps 
in a given mass range, two types of contributions exist. Fol- 
low mass clumps, which arc numerous, the variance associ- 
ated with the flux is expected to be small (i.e. cr c \/ <p c \ <C 1). 
In this case, we can spare the effort of averaging many con- 
figurations and directly compute the flux from a single real- 
isation. Conversely, as the mass of the sub-halos increases, 
the associated number of clumps decreases, so that the vari- 
ance finally become sizable. A threshold mass M t h needs 
to be specified, below which the contribution to the total 
variance ct'°* can be neglected: only sub-halos that have 
masses M c \ > M t h need to be calculated for all samplings. 
The value of M t h is discussed in App. 

For one sample, the total annihilation flux observed in 
the solar neighbourhood may be rewritten as 



6 tot 



'high, 
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where the quantities 



'low 



Mi 



— S I d 3 x / dM c i 0(x x )n(M c i, r) 



"high 



Psub 
'sub \ PO 

Mi>M th 



2 



(x')d 3 x' 



(73) 
(74) 



are the contributions from the low-mass sub-halos compo- 
nent and the high-mass sub-halos component, respectively. 
The number density of clumps was defined in Eq. and 
the flux from a single clump <j>i is given by Eq. (|51 
When taking into account all realisations: 



'high/ 



where (.)„ denotes the average over n 3> 1 realisations of 
the spatial distribution. This leads to: 




and 



cl } = "Alow + (0high) 



-_tot _ 
O c \ _ <Thigh 



6 tot } 



6 to t) > 



(75) 



(76) 



where (Thigh is the variance associated to high-mass clumps 
(<ti ow is neglected as underlined above). 

The total flux and variance are now given by 



and 



= 9am + 


01ow + 


Ctot 


Chigh 


(0tot) 


('/'tot) 



'high/ 



(77) 



(78) 



Thanks to this reasoning, a simple picture emerges, and 
a qualitative behaviour of the expected variance nicely com- 
plements the discussion from the previous method. The 
highest-mass sub-halos, which are rare, carry all the vari- 
ance of the total flux. In App. IA.31 they are found to be 
in the mass range M c \ > 1O 7 M0. Thus, as soon as the 
integrated luminosity of lower-mass clumps (that depends 
on a m ) is much larger than that of the high-mass ones, 
the variance of the total clump contribution is expected to 
vanish. According to Fig. [5J such a situation will occur for 
a m si 1-9. Actually, even for a m < 1.9, the variance will 
still be significantly decreased, because, as we already un- 
derlined, the local smooth contribution dominates over the 
clump one (see Fig. [5]). 

6. Results and discussion 

Fluxes, boost factors and associated variances have been 
calculated for both positrons and anti-protons using a semi- 
analytical approach (Scct. lBT^)) and, for the sake of compar- 
ison, MC simulations (Sect. [573]) . 

The coming results are based on a fiducial model 
for the injection of antimatter in the Galaxy, which al- 
lows a WIMP-model-independent analysis. In practise, for 
positrons, a monochromatic line of 200 GeV is injected at 
a rate assumed to be proportional to the squared density 
of DM in sources. In order to recover realistic orders of 
magnitude, especially for fluxes, we will also suppose that 



those positrons originate from, e.g. not-s-wave-suppressed 
annihilations of WIMPs at rest, with masses of 200 GeV 
and annihilation cross-section (av)/2 = 3 x 10 -26 cm 3 .s _1 
(for instance, Dirac fermions/anti-fermions with only trilin- 
ear couplings to e~0 + /e + </V, where <j) + l~~ would be some 
exotic — conjugate — charged scalar fields). Besides, because 
anti-protons do not loose energy, we took their injection 
spectrum to be constant dN/dT — lGeV^ 1 between kinetic 
energies 0.1-200 GeV (any spectrum could have been taken, 
as it can be factorised out). One can easily guess what the 
results for any injection spectrum would be (originating 
from hadronisation or fragmentation processes for instance) 
by a mere rescaling. In this case, the WIMP properties can 
be almost the same as for positrons: fluxes have been com- 
puted using a Majorana WIMP with a mass of 200 GeV, 
and an annihilation cross-section of (av) = 3 x 10 _26 cm 3 /s. 

Before going into the details of the studied configura- 
tions, we show in Fig. [5] the extreme cases that we obtained 
for both species (with the medium set of propagation pa- 
rameters). The first line panels are plots of the smooth and 
sub- halo fluxes and the resulting effective boosts, with asso- 
ciated 1 — (7 statistical contours. The second line panels are 
the same plots, but for anti-protons. The maximal config- 
uration is given by: largest a m (2), cuspiest sub- halo inner 
profile (Moore), smallest M m ; n (10~ 6 M Q ), spatial distribu- 
tion according to the smooth NFW profile, and the B01 con- 
centration model. The minimal configuration is the reverse: 
smallest a m (1.8), flattest inner profile (NFW), greatest 
M m i n (10 6 Mq), and smallest local number density (cored 
isothermal profile). The intermediate is close to the refer- 
ence configuration, given in Sect. [3~l~2l and takes the most 
likely values of parameters according to N-body simula- 
tions (except for Af m i n , of which the used reference value is 
10~ 6 AfQ; and for the spatial distribution of sub-halos which 
tracks the smooth NFW profile). 

From this figure, we see that the boost factors obtained 
are functions of the energy and lie between 1 and 20, with 
small statistical uncertainties. Such a range has to be taken 
as that of theoretical uncertainties affecting the DM distri- 
bution in the Galaxy. From the approximate Eq. (|68p . the 
asymptotic values obtained are also 1 and 20 (neglecting 
the density fraction /©), which arc in excellent agreement 
with our full results. Before going into deeper details, it is 
worth emphasising that the maximal value of ~ 20 is as 
large as unlikely, as already discussed in Sect. [2] For com- 
pleteness, we have checked our results with MC simulations 
(see Sect. I5.3[) . In Fig. we show that the agreement be- 
tween the MC and the semi-analytic calculation is excellent 
for anti-protons, up to a few percen ts. It is the same fo r 
positrons, as alrea dy demonstrated in lLavalle et al.l (|2007t ) ; 
iBrun et all (|2007h . 

Remembering that the whole sub-halo flux reads </>*° t cx 
(V c i(£)m) x ((G) v), it makes sense to gather the impact of 
the various ingredients into two main physical classes. 

Mass-related effects (N c \(£)m): these are encoded in the 
mean value and the variance of £, as defined in Eqs. (|6"2"]l 
and (|65[) . The relevant parameters to discuss (see 
Table ^ are the minimal mass M m - m of sub-halos, the 
logarithmic slope of the mass distribution a m , the mass- 
concentration model and the inner profile. 

Space-related effects ((G) v)'- these are encoded in the 
mean value and the variance of G as defined in Eqs. (|6ip 
and (|64|) . which depend on both the propagation model 
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Fig. 6. Extreme cases for the DM configurations: sub-halo antimatter fluxes associated with the maximal, intermediate and minimal 
DM configurations (medium set of propagation parameters). Left/right: fluxes/boosts and corresponding 1 — a contours. Top/bottom: 
positrons/anti-protons. See the details in the text. 



(through the propagator and the propagation parame- 
ters) and the spatial distribution of sub-halos, besides 
of course the antimatter species. 



More details on this classification can be found in the ap- 
pendix (see App. [Bj. We first discuss the mass-related ef- 
fects. We then focus on each antimatter species, for which 
we explain space-related effects, before describing the whole 
consequences of dumpiness on fluxes and boost factors. 
When discussing positrons, we will comment on the so- 
called HEAT excess. 



6.1. Mass-related effects 

Once the space distribution of sub-halos and the propaga- 
tion model are fixed, the propagator mean value (G)v of the 
sub-halo flux Eq. (|STI)) is fully determined, as well as its sta- 
tistical fluctuation <jq. Hence, provided the WIMP model 
is also fixed, the only differences from one sub-halo config- 
uration to another will be the averaged total amount of an- 
timatter yielded by clumps, given by the integrated clump 
luminosity L c \ = N c \ x (£) a/, and its associated fluctuations. 
Such a quantity depends on two parameters only: M m i n and 
a m (plus the concentration-mass relation, plus the choice 
of the inner sub-halo profile) . A decrease of M m i n enhances 
the total number of clumps in the Galaxy, and an increase 
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Fig. 7. Ratio of MC over semi-analytic results for the anti- 
proton flux and associated variance. 



of a m raises the relative density of light compared to heavy 
objects (and the total luminosity accordingly because the 
clump number density is normalised with respect to the 
heaviest clumps, as given in Eq. [17|) . More precisely, we 
find the luminosity to approximately scale with A/ m j n only 
logarithmically (see the details in App. IB.1[) like: 

L cl cx M*?" 1 x ln(M max /M min ). (79) 

Therefore, we do not expect a large variation when span- 
ning (a m ,Af m i n ) from the minimal (1.8,1O 6 M0) to the max- 
imal (2,10 _6 M Q ) parameter sets. Actually, we find the total 
luminosity to drop by a factor of ~ 40 only, while running 
the number of clumps over 13 orders of magnitude in the 
meantime (the above approximate expression sightly over- 
estimates the difference, leading to a factor of ~ 150). For 
completeness, an additional factor of 10 appears if one also 
considers inner Moore profiles for clumps (see Sect. I3.3TT) . 
whereas a factor of 1/10 results from using an ENS01 mass- 
concentration relation. 

Regarding the pure mass-related relative fluctuations, 
given by N C \~ 1 ^ 2 ((T^/ (£)) (see Eq. lM)) . we would naively ex- 
pect them to significantly deplete when decreasing M m ; n , 
which enhances the total number of clumps. Furthermore, 
since the relative luminosity of light clumps is raised by in- 
creasing a m , the relative variance should be significantly 
reduced accordingly. Nevertheless, interestingly, we find 
the mass-related relative fluctuations to roughly scale like 
(ln(M max /M m i n )) _1 , and to vary only in the range 0.1- 
10% when spanning over (a m ,M m ; n ) from minimal to maxi- 
mal parameter sets (see above). The physical interpretation 
is the following: as M m j n goes up, the number of clumps 
decreases accordingly, but in the meantime, the intrinsic 
clump luminosity (oc £), which is fixed by the mass, fluc- 
tuates much less from clump to clump; there is a trade-off 
between shrinking the statistical sample and reducing the 
phase space, so that the relative mass-related variance re- 
mains almost constant. Taking a Moore inner profile does 
not affect the relative mass-related variance, while adopt- 
ing the ENS01 mass-concentration relation increases it by 
a few (see the third line panels of Fig. [5]). 

As a consequence, the global flux relative variance given 
in Eqs. (|63p and (|66|) . should vary over two orders of magni- 
tude at most, once the GCR propagation is fixed and at a 



given energy. Indeed, as mass-related uncertainties are al- 
most always greater than space- related ones (see App. IB.lj) . 
the dominant contribution is the space-mass crossing term 
N c i~ 1 / 2 (a£/(£ t ))(o-g/(G}) 1 so that the global relative vari- 
ance encompasses values in the range (0.1 — 10%) x (cg/(G)) 
when varying a m and M m i n from extreme configurations. 

From the previous statements together with the lumi- 
nosity profiles already discussed and shown in Fig. [5l scan- 
ning over the most likely mass-related parameters is un- 
likely to make the sub-halo contribution dominate over the 
smooth flux, except for extreme configurations combining 
the B01 concentration model, Moore inner profiles, large 
a m and very small M mnl . 



6.2. Positrons 



6.2.1. Space-related effects for positrons 

The space-related effects for positrons come through the av- 
eraging of the propagator (G e+ )v over the sub-halo spatial 
distribution. We summarise here a more detailed discussion 
that will be found in the appendix (see App. lB.2"|) . The rele- 
vant scale is the propagation scale Xd that depends on both 
diffusion and energy loss processes for positrons. Xd is ob- 
viously larger for larger diffusion coefficients, and smaller 
when the detected energy gets closer to the injected energy. 
Since it is of the order of kpc, we can safely focus on local 
quantities. Actually, (Q e )y encodes an effective detection 
volume bound by Ad and weighted by the clump spatial 
probability function dV/dV(r) in the solar neighbourhood. 
In the limit of infinite 3D diffusion, and when the prop- 
agation length is small enough, we find in App. IB. 21 that 
(G e+ ) ~ (r £ /e 2 ) x dV{R®)/dV . Hence, the averaged prop- 
agator increases linearly with the local value of the clump 
spatial probability function. As dV / dV {Rq)(p & / M^) for 
the reference case (clumps are spatially distributed accord- 
ing to a cored isothermal profile), we see that given mass- 
related parameters, a configuration in which the clumps 
track the smooth profile will give a higher flux. 

Regarding the pure space-related relative variance for 
a single object <Jg/(G), we find it to scale like (X D x 

dV{RQ)/dV)~ 1 / 2 , thus, decreasing when the effective de- 
tection volume or the clump local spatial probability in- 
crease (detected energies much lower than injected ones). 
When taking the whole contribution, an additional factor 

1/9 

of N c \ 1 reduces the global variance, and the picture be- 
comes very simple: the relative space-related variance scales 

— 1/2 

like N b , one over the square root of the number of 
clumps contributing to the signal at the Earth. It is maxi- 
mal at high energy for positrons. 

To summarise, the space-related contribution for 
positrons increases with the diffusion coefficient, and with 
the clump local space probability function. The relative 
space-related variance decreases when the propagation 
length raises (at low energy for positrons) , because a larger 
number of sub-halos can contribute to the signal at the 
Earth. 
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6.2.2. Overall effect on the positron flux: boost factor 
estimate 

Taking this fiducial injection model, we assess the differ- 
ent effects and draw four typical plots, which will compose 
four specific panels in the next figures, from left to right: 
positron flux, relative flux variance, boost factor, relative 
boost variance (as functions of the positron detected en- 
ergy). 

Figure [5] illustrates the mass-like effects, whereas FigJ5] 
show the space-like ones. 

— Mass-related effects 

In the first line of Fig. [51 we vary the minimal mass of 
the sub-halos, in other words the cut-off of the mass distri- 
bution. We actually compare three configurations by tak- 
ing M m i n = 10 -6 , 1 and 10 6 M@, the remaining param- 
eters being those of the reference configuration given in 
Sect. 13.11 The top left panel shows the whole contribution 
of sub-halos to the positron flux with the associated 1 — a 
contour, as well as the smooth contribution. We see that 
varying the minimal mass mainly influences the variance, 
while the mean values predicted for the flux remain of the 
same order of magnitude (only a factor of ~ 5 in flux be- 
tween M min = 1CT 6 and 1O 6 M ). The flux ratios of the 
three configurations are plotted in the second panel of the 
first line, taking the reference 0_ 6 = <f>(M m i a = 10~ 6 M Q ). 
This is due to the fact that the product N c i x (£) is almost 
independent of the minimal mass in this case. Should a m 
have taken a value different than 1.9, the mean contribu- 
tion of clumps would have a much stronger dependence on 
M m ; n . Nevertheless, the relative flux variance is different 
between the three configurations, as also shown in the bot- 
tom left panel. Actually, this comes from the total number 
of sub-halos, which is strongly depleted when M m ; n is in- 
creased (oc M m i n ~ am ). Therefore, the statistical flux vari- 
ance, which scales like 1/^N C \, is increased accordingly. 
Another effect comes from the energy loss of positrons. 
While the detected energy gets closer to the injected energy, 
the diffusion volume decreases, as explained in Sect. 16.2.11 
and the number of clumps effectively contributing damps 
in the same way. Hence, the relative variance is enhanced 
when getting closer to the injected energy. However, the 
whole sub-halo contribution is finally far below the smooth 
flux, by about two orders of magnitude. This translates to 
an effective boost of B e ff(Ed) — (1 — /©) 2 ~ 1, with a very 
small variance, because it is also diluted by the smooth 
component. 

In the second line panels of Fig. [5J a m is varied, giv- 
ing three different mass configurations: 1.8, 1.9 and 2.0. As 
expected, the flux due to sub-halos is affected, and pre- 
dictions slightly spread within one order of magnitude (a 
factor of ~ 30 between 1.8 and 2.0). The relative flux vari- 
ance is lower for large values of expected, because 
this increases the total number of sub-halos, more precisely 
the lighter ones. Nevertheless, varying a m within the refer- 
ence configuration is not enough for sub-halos to strongly 
dominate over the smooth contribution: the averaged boost 
factors associated with the three examples lie around unity, 
even when getting closer to the injected energy (the max- 
imum value is 1.08 for a m = 2.0), with small statistical 
uncertainties. 

The third line of Fig. [8J shows the consequences of vary- 
ing the mass-concentration relation and the inner sub-halo 



profile. The reference model, which is inner NFW + B01, is 
compared with NFW+ENS01 (less concentrated sub-halos) 
and with Moore+B01 (more cuspy sub-halos). As expected, 
the flux obtained with the ENS01 concentration model is 
far below the reference one, by a factor of ~ 20, whereas 
the Moore sub-halos gives ten times more signal (this can 
also be seen from Fig. [3]). These ratios are constant with the 
detected energy, as they are characterised by the ratios of 
(£)'s. Again, we see that the expected boost factor is again 
negligible in all cases, around unity (~ 1.13 for the best 
case, i.e. Moore inner profile + B01). Nevertheless, the in- 
crease of the variance associated with the latter happens at 
lower energies than previously, because the probability for 
a single clump to contribute more than the smooth compo- 
nent becomes sizable at farther distances. 

— Space-related effects 

The first line panels of Fig. [5] illustrate how GCR propa- 
gation can strongly influence the predictions by using the 
three models of Table [4] with the reference DM configu- 
ration. Differences grow when the detected energy is far 
below the injected one. This is mainly due to the change in 
the thickness of the diffusive slab. At energies close to the 
injected energy, the volume probed is very small, so prop- 
agation is not sensitive to the slab boundaries anymore. 
We see that for positrons, there is no huge differences be- 
tween the maximal and medium propagation models, while 
the minimal one strongly depletes the positron flux at low 
energies. Indeed, the characteristic propagation length for 
positrons (a few kpc) is almost always contained in the 
maximal and medium slabs, whereas it is not the case for 
the minimal one. Anyhow, even the maximal propagation 
set is not enough to boost the sub-halo positron flux above 
the smooth contribution. Indeed, the values obtained (see 
the left panels) stack to unity, again with a small statistical 
uncertainty. 

Finally, the second line panels of Fig. [5] show the effect 
of changing the spatial distribution of clumps, from the ref- 
erence cored isothermal to a situation in which they track 
the smooth NFW profile. For completeness, we do the ex- 
ercise for both inner NFW and Moore profiles. In the left 
panels, we see that the sub-halo flux is enhanced when they 
track the smooth profile, of about one order of magnitude 
in this case. The effect is obviously stronger when an inner 
Moore density is taken, for which another order of magni- 
tude arises. Nevertheless, the boost factors do not obey the 
same hierarchy. This is due to the way the smooth com- 
ponent is normalised when clumps are added. Indeed, we 
chose to readjust the smooth density by a factor (1 — /©), 
where the fraction density /© is defined in Eq. (|2"2")h in or- 
der to get a constant local density p Q . When clumps track 
the smooth component, their number density is enhanced 
in the local neighbourhood, so that /© is enhanced accord- 
ingly. As the boost factor remains around (1 — /©) 2 (the 
sub- halo contribution is negligible), except when consider- 
ing inner Moore profiles, the case for which clumps track 
the smooth component is worse. However, taking an inner 
Moore profile for clumps gives a higher flux than the smooth 
alone, and the mean boost factor can reach an asymptotic 
value of ~ 2 — 3. 

In summary, we have shown, by extensively playing 
with the maximum number of available parameters, that 
sizable boost factors to the positron flux are unlikely to 
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Fig. 8. Mass-related effects on the positron flux (first column) with a focus on the corresponding relative uncertainty (second 
column), for positron lines of 200 GeV injected at sources, with a rate corresponding to standard values of WIMP annihilation, 
and the associated boost factor with its relative variance (third and fourth columns). The plain contours account for one standard 
deviation. First row: Effect of changing the minimal mass of sub-halo s, M m i n = 10~ 6 , 1 and 10 -6 Mq. Second row: Effect of 
changing the logarithmic slope a m of the sub-halo mass function, with a m = 1-8, 1.9 and 2.0. Third row: Effect of changing the 
sub-halo inner properties ; the inner profile is taken to be either NFW or Moore, and the concentration model varies from B01 to 
EN SOI. 



arise from dumpiness. There could be situations in which 
a single sub-halo would be close enough to the Earth to 
dominate over the smooth component, but within the most 
reasonable modellings, the probability for this to happen is 
vanishingly small. Nevertheless, in order to provide more 
optimistic scenarios, even if less realistic according to the 
standard values of the parameters, a model characterised 
by 10 -6 Mo sub-halos with inner Moore profiles, with con- 
centrations described by the B01 relation, spatially tracking 
the smooth DM density, would yield a mean boost factor 
whose asymptotic value would be around 3 for a logarith- 
mic slope a m = 1.9 (see lower panels of Fig. [5]). Taking a m = 
2 leads to a boost of ~ 20 (see Fig. [5]). This is the most op- 
timistic estimate that we can provide so far, but also the 
most unrealistic. Note finally that although the primary 
fluxes may vary by 2 orders of magnitude due to uncer- 
tainties in the propagation parameters (see Table [4] and 
upper panels of Fig. [5]) , the resulting average boost factors 
are unaffected because they are defined as flux ratios; the 
variance is nevertheless larger when the GCR horizon is 
reduced (min configuration). 



6.2.3. Comments on the positron excess 

The HEAT experiment results for the 1994 flight hinted 
at th e possibility of an excess of positrons near 8 
GeV (jBarwick et al.|[l997|) . which could not be explained 
by a p urely secondary produ ction mechanism (|Coutu et al.l 
ll999D . lBaltz fc Edsioi(|1999D then found that neutralino an- 
nihilation could account for the missing flux providing that 
boost factors are larger than six; at that time, these authors 
estimate realistic boost factors to fall in the range B e g < 
100 — 1000. Note that such high boost factors would be 
ruled out in the present study. However, la ter on, combining 
both 1 994 and 1995 HEAT balloon flights. iDuVernois et al.l 
()200lh concluded that the positrons flux was consistent 
with a secondary origin. Results from the M ASS91 balloon- 
born e magnetic spectrometer above 7 GeV ([Grimani et al.l 
|2002| ) do not provide a definitive answer either. As empha- 
sised by these authors, very high energy ~ 100 GeV mea- 
surements arc probably necessary to positively conclude for 
standard or exotic mechanisms. Finally, from the most re- 
cent data coming from the HEAT 2000 flight. iBeattv et al.l 
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Fig. 9. Space-related effects on positrons, with the same panel organisation as in Fig. First row: Effects of varying the GCR 
propagation modelling, by using the three propagation sets of parameters of Table rj[ Second row: Effect of varying the space 
distribution of sub-halos, going from the cored isothermal space distribution to a case in which sub-halos track the smooth NFW 
component (for completeness, we also take an example where we also modify the inner sub-halo profile by taking a Moore instead 
of a NFW). 



contribution 



(|2004f ) cautiously conclude that a primary 
above a few GeV can still not be ruled out. 

Given these observations, several subsequent studies 
have focused on finding a good DM candidate to ex- 
plain this possible excess. We do not wish to com- 
ment here on the best candidate, but rather survey 
the boost factors used in the studies. For example, for 
SUSY candidates, boost f actors of 2.7 and 3.9 were 
used in iKane et alj (|2002bD . values in the range 30-100 
2002T). 



m 
in 



Baltz et al 



from 



m 

small to large 



boost factors 



.Kane et al.l (|2002ah . in t he range 1-5 in Irlisano et alj 
(120061). and around 1 00 in ICumberbatch fc Silkl (|2007l) . 



iMambrini et alj (|2006f ) favoured boosts of 5-10 to accom- 
modate the expected measurements of PAMELA, for SUSY 
models w ith non-universal scalar and gaugino masses. For 
KK DM, iHooper fc Kribsl (|2004) found a range of 10-30. 
The boost factor used to fit the data depends of course on 
the WIMP candidate considered and its mass. 

It appears that most of the models found so far to 
match the positron data require mild to significant boost 
factors. Such boost factors are disfavoured by our results if 
the clump parameters fall in the large ranges taken in this 
study. A high energy feature in the positrons data could 
still be an important clue to DM indirect detection, and 
it would be interesting in forthcoming studies to scan, e.g. 
the SUSY parameter space looking for mod els matching 
the da ta without boost factors. According to lAsano et alj 
(|2007| ). the little Higgs model provides good options for de- 
tectability by the AMS-02 experiment, but could be short 
for PAMELA. In addition, it is worth noting that if some 
X 2 -likc searches for clump signatures are performed in the 
coming positron data, it will be very important to take the 
energy dependence of any boost factor into account, as soon 
as it is invoked. 



To conclude, although having no boost factors may be 
less interesting for SUSY theories to explain the data, any 
result that will be obtained when comparing to forthcom- 
ing data, if an excess is confirmed, will be more robust 
if no boost factor (an additional unknown parameter till 
now) is invoked. We recall that the two main uncertainties 
for WIMP annihilation induced antimatter signals are the 
propagation parameters in the Galaxy (a factor < 100) and 
the local DM density (a factor < 2, that shifts to 4 in terms 
of annihilation rate). 

6.3. Anti-protons 

First of all, it is worth quoting that contrary to positrons, 
for which an excess is still not understood, anti-proton 
present measurements are now well ac counted for by purel y 
standard secondary production fe.g. iDonato et aT1l200l[ ). 
This means that there is no need of DM, and obviously 
of any clump to fit the data: the present data has to be 
consider as an upper limit for the DM contribution. Things 
could change with the future results of PAMELA and AMS- 
02 at higher energies. 

The flux enhancement for anti-protons has fea- 
tures different from posi trons , as already stre ssed by 
iBringmann fc Salatil (|2007h and lBrun et~aT1 (|2007ft . This is 
mainly due to propagation, which is quite different from 
the positron case. Indeed, anti-protons do not lose energy, 
and can experience spallation processes and wind convec- 
tion along their travel to the Earth, which occurs to be 
dominant at low energy Nevertheless, as for positrons, the 
same classes of physical effects can be discussed. 

Regarding the space-related effects, the comments are 
the same as for positrons but with a reversal energy point of 
view (see App. |B~3| . The relative variance is then maximal 
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at low anti-proton energy. Note, however, that the three 
sets of propagation parameters give separate absolute fluxes 
(decreasingly according to max, med and min), whereas 
max and med configurations give about the same fluxes for 
positrons. 

6.3.1. Overall effect on the anti-proton flux: boost factor 
estimate 

We now discuss the origin of systematic differences when 
varying the DM configuration as well as the propagation 
modelling. Figure [TU] illustrates effects that are of mass 
type, while Fig. [TT] shows the space-related ones. They are 
presented the same way as for positrons. 

For the mass-type category, the comments are exactly 
the same as those for positrons, and are already discussed 
in Sect. 1X1 

Regarding space-related consequences, the picture is the 
reversal from that of positrons, and the conclusion are the 
same as for positron, given the energy axis is read inversely. 

In summary, the maximum boost factor occurs at low 
anti-proton energies when clumps are spatially distributed 
according to the smooth profile, and when they have an 
inner Moore profile. But even in this (disfavoured) config- 
uration, the asymptotic mean value of the enhancement 
factor is < 3. 



7. Summary and conclusions 

dumpiness is a robust prediction of hierarchical structure 
formation, where the bottom-up growth of structures is a 
consequence of the ACDM cosmology. Many issues remain 
about clumps, beside their survival, such as their number 
density, their mass and spatial distribution, as well as their 
intrinsic properties. In recent years, high resolution N-body 
simulations succeeded in tracing their gross features, al- 
though we arc still far from a definite answer for several key 
parameters. Indeed, whereas numerical experiments now 
converg e to a level better th an 10% over wide dynamic 
ranges (|Heitmann et al.ll2007f ). the inclusion of baryons in 
hydrodynamical simulations, which could strongly modify 
the inner properties of sub-halos, remains a very difficult 
task. 

The impact of the DM inhomogeneitie s on DM annihi- 
lation in the Galaxy was first underlined in lSilk fc Stebbinsl 
(|1993| ). With the achievement of cosmological N-body sim- 
ulations in the last decade, and due to their positive results 
on dumpiness, many papers have subsequently focused on 
the consequences for the phenomenology of SUSY indirect 
detection in sev eral channels, involving 7-rays and antimat- 
ter GCRs (e.g. iBergstrom etaTl 119991) . The effects for 7- 
rays have been extensively studied, whereas diffusion pro- 
cesses make the study much m ore difficult for antimatter 
GCRs. lBerezinskv et al.l (|2003l ) carried out a detailed study 
of boost factors for 7-ray signal (so the conclusions also hold 
for neutrinos), and concluded to a factor of 2 to 5 for dif- 
ferent density profiles in a clump. However, whether or not 
primary antimatter fluxes may be boosted by dumpiness 
is an important issue for the interpretation of forthcoming 
data, especially for positrons (see the discussion below), 
and for putting more robust constraints on new physics. 

In this paper, we have inspected the boost factors 
for antimatter GCRs as deeply as possible and the re- 
sults are mostly independent from any WIMP model, in 



the context of cosmological substructures. To this aim, 
we have followed th e semi-analytical method proposed in 
lLavalle et all (12007), a l ready used by iBringmann fc Salatil 
( 20071 ) and lBrun et al.1 (|2007[ ) to compute the boost factors 
for positrons and anti-protons, but in a scenario in which 
DM inhomogeneities are due to the presence of intermedi- 
ate mass black holes. We have also cross-checked our results 
with a more time-consuming MC simulation. 

A full DM model was defined with (i) the host halo 
smooth DM profile, (ii) inner properties, minimal mass, 
space and mass distributions of sub-halos, and (iii) the 
mass-concentration relation. The reference set was cho- 
sen accordingly, with (i) NFW (kept fixed throughout the 
paper, see Table [J), (ii) inner-NFW + 10~ 6 M Q + cored 
isothermal + a m = 1.9, and (iii) B01 concentration. We 
then extensively spanned over some extreme values of each 
parameter. We normalised all configurations so that the av- 
erage local DM density was kept constant p t ot(-R©) = Pq — 
0.3 GeV cm -3 . For completeness, we also considered three 
sets of GCR propagation modelling, mainly characterised 
by the size of the diffusive halo, and by the amplitudes 
of diffusion and convection processes; all of them decrease 
from max to min settings (med being the reference). 

We found our results to depend on two main classes of 
effects, namely mass- and space-related. The mass-type ef- 
fects characterise the full amount of antimatter produced in 
clumps through WIMP annihilations, given by the product 
N c i(£)m (see Sect. I6.1[) . and its variance; they depend on 
the number of clumps (fixed by M m - m at a given a m ), their 
inner profiles (NFW or Moore) and mass distribution (a m ). 
The space-related effects describe the average probability 
for the produced antimatter to reach the Earth, and is de- 
termined by (Q), and the associated variance; they therefore 
encode the whole spatial information, i.e. the propagation 
averaged on the clump space distribution, as well as the en- 
ergy dependency. Our results are the following, some trends 
being consistent with some found in previous studies: 

— larger a m values lead to larger luminosity of clumps 
(hereafter C c i oc N c i(£)m), because N c i increases faster 
than (£)m decreases; 

— diminishing M m j n values lead to larger C c \ if a m > 1.9, 
otherwise this has no effect on the average luminosity 
(massive clumps dominate the signal in average); 

— the luminosity of a clump of a given mass with an inner 
Moore profile is 10 times that obtained for the same 
clump with a NFW profile (which can be analytically 
handled) ; hence, the whole clumpy flux can be evaluated 
from a NFW profile and rescaled to any profile very 
easily, being for instance 10 times larger for a Moore 
profile; this is independent of the concentration relation 
assumed and the species considered; 

— the maximum flux is obvioulsy obtained for maximal 
concentrations, here the B01 c v i r — Afvir relation, and 
is lower for the ENS01 case (respectively a factor of 
~ 3 and ~ 30 lower — mass dependent — for M c \ = 
10 10 Mq and M cl = 10~ 6 M ); 

— spanning extreme ranges of parameters leads to the 
same conclusion that no sizable enhancement is ex- 
pected from dumpiness in average, asymptotic values 
of B c ft varying from (1 — /q) 2 ~ 1 to ~ 2.5; 

— boost factors for GCRs depend on energy, which char- 
acterises the size of the effective volume of sensitivity 
that can be probed around the Earth; 
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Fig. 10. Mass-related effects for anti-protons ; a toy flat spectrum is injected at the sources, at a rate corresponding to standard 
values of WIMP annihilation — same panel organisation as for positrons (see Fig. \Ej). First row: Varying M m i . Second row: 
Varying a m . Third row: Varying the sub-halo inner properties. 



— the relative statistical uncertainties on the whole sub- 
halo fluxes and on the corresponding boost factors are 
dominated by mass-related effects — linked to the in- 
ternal properties of clumps and to their mass distribu- 
tion — as soon as small clump masses are considered 
(< IMq), and are hugely enhanced due to the mixing 
with the (weaker) space-related ones — linked to the 
spatial distribution of clumps and to propagation - 
especially when the GCR propagation scale gets slen- 
der and slender, i.e. for low energy anti-protons, and 
detected (high) energies close to the injected one for 
positrons (while the absolute fluxes do depend strongly 
on propagation); 

— statistical uncertainties are small at low energies for 
positrons and at high energies for anti-protons (when 
the propagation scale is larger) . 

As a first and mere conclusion concerning the whole 
sub-halo flux, the combination of smooth-like space dis- 
tribution, small M m i n , large a m , large concentrations, and 
very cuspy inner profiles will obviously lead to the largest 
mean flux (in the statistical sense of averaging over many 
DM outcomes for a given configuration). The maximum 
set of parameters that wc considered for the sub-halos 
was defined by a spatial smooth-like NFW distribution, 



-Mmin = 10~ 6 M Q , a m = 2, a Moore inner profile (cx J" -1 ' 5 ), 
and the B01 concentration model. This led to asymptotic 
values of boost factors around 20 with small statistical er- 
rors (see Fig. [6]). Conversely, a cored isothermal spatial dis- 
tribution, large A-/ m j n , small a m , small concentrations (here 
ENS01) and mildly cuspy (e.g. NFW) inner profiles will 
lead to the smallest mean flux. We remind the reader that 
a very simple and straightforward way to estimate whether 
or not sub-halos may enhance the DM contribution to the 
antimatter fluxes is to verify whether the condition given 
by Eq. (JSHD is fulfilled. 

Furthermore, note that the only relevant parameters for 
estimating the variance on the sub-halo flux are the mini- 
mal mass of clumps, their mass distribution and their local 
number density, beside the propagation length of the GCR 
species that defines an effective detection volume. 

These re s ults are in agreement with those of 
iHooper et alJ (j2004h . who concluded, though mostly qual- 
itatively, that it is very unlikely that significant boost fac- 
tors occur for positrons. Our results, however, are more 
quantitative and detailed, apply to both anti-protons and 
positrons for any set of propagation parameters, and en- 
compass the single configuration (NFW-ENS01) used in the 
above paper for clumps. This is also consistent with the re- 
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Fig. 11. Same than Fig. but focusing on space-related effects for anti-protons. First row: Varying the propagation modelling. 
Second row: Varying the space distribution (and the sub-halo inner profile). 



suits of iBerezinskv et ail (|2003|) for 7-rays, where slightly 
larger boost factors were found: such a difference is ex- 
pected because, as already underlined, 7-rays arc integrated 
along the line of sight (we remind that the luminosity of 
clumps dominates over the smooth distribution one beyond 
a few tens of kpc from the Galactic centre) , instead of in- 
tegrated inside a more local volume like for GCRs. 

It is interesting to ask the dependence of our calculation 
on the mass resolution achieved so far in N-body simula- 
tions, which we have referred to in defining our param- 
eter sets. Indeed, the N-body numerical results arc only 
valid at the spatial scale associated with the test parti- 
cle mass, so that the Vlasov limit may not be reached 
at the smallest scales considered here (see discussion in 
Sect. [2]). Therefore, extrapolations of the physical proper- 
ties of sub-halos down to 1O _6 M should always be taken 
cautiously, even if some numerical studies were able to sur- 
vey su ch small systems at high rcdshifts (jDiemand et alJ 
l2005al h Besides, even if DM sub-halos of 1O _6 M wander in 
the Galaxy, we can actually not know anything about their 
characteristics. Nevertheless, theoretical arguments based 
on the (inflation-motivated) scale invariance of the DM 
power spectrum down to the free streaming scale set by par- 
ticle physics, the theoretical understanding of the DM mass 
function, the current knowledge, would it be far from com- 
plete, of hierarchical structure formation, and some numeri- 
cal studies on the survival of very small sub-halos, somehow 
guarantee that our choice of parameter ranges is rather rea- 
sonable and sufficiently large to encompass a wide field of 
possibilities. Therefore, while this strongly asks for more 
detailed studies of the smallest DM structures, our results 
should also be taken as general statements that describe the 
effect of each considered parameter on boost factor predic- 
tions. 

One could recover a sizable (energy dependent) boost 
factor by considering either a sub-halo which would be 



very massive (> 10 7 Af Q ) as well as very close to the Earth 
(< 1 kpc, see e.g. Sect. I6.2.ip . or very cuspy inner pro- 
files combined with a significant local abundance of sub- 
halos. The latter case would correspond to a clump con- 
figuration given by a smooth-tracking spatial distribution, 
M min < 10~ 6 Mq, a m > 2, concentration > B01, and r - ' 9 
inner profiles with /3 > 1.5. Such a situation is very improb- 
able given the current theoretical results of gravitational 
collapse or mass function studies, and also considering the 
most likely configurations of dumpiness found among N- 
body results. Thus, it may not be taken as a natural predic- 
tion of structure formation. Regarding the former case, we 
stress that the statistical probability to find such a massive 
object in the solar neighbourhood is vanishingly small (such 
masses are now well resolved in numerical simulations, and 
are not expected to be numerous). Moreover, some obser- 
vational constraints might exist on the presence of such a 
massive and close object. Anyway, the calculation of fluxes 
originating from a single nearby source is straightforward, 
and one can very easily model its required features. The 
price to pay would be to invoke some kind of Galactic lot- 
tery in order to explain why a single clump would wander 
here and now. 
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Appendix A: Effective propagation volumes and 
the threshold mass M t h 

On the one hand, the spatial distribution and the num- 
ber of clumps in each mass decade is known. On the other 
han d, the propagation properties define effec tive volumes 
V eS (|Taillet fc Maurinll2003l ; fMaurin fc TailletJl2003D . which 
enclose and pre-select the sources contributing to the flux. 
This is all we need for a quick estimate of the variance on 
the fluxes (for the clumps in that given mass range), hence 
the estimate of the threshold mass M t h- Effective volumes 
V e g are reminded in App. IA.fi and Mth is given in App. IA.2l 
In App. IA.3l we show that M t h can be set to a higher value 
than that derived from the effective volumes. 

The method is general, and is discussed below for a m = 
1.9 (Sect. [2~3Jl Eg. [14)) and the cored distribution of clumps 
(Sect.EXa Eq.Mj. 



A.l. Effective volumes V c k 

All DM sources beyond the boundary z = L (size 
of the diffusive halo of th e Galaxy) can be safely dis- 
carded (jBarrau et al.l l2002f ). Furthermore, in a diffusive 
process, a source located at a radial distance r gives a 
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negligible contribution if 77 = r/L is larger than a few 
(jTaillet k Maurinll2003l ; IMaurin fe Tailletll2003l ) . These two 
boundaries generate a cylinder V e s = irrj 2 L 2 x 2L; sources 
out of it may be considered to add negligible contribution 
to the total fkuQ 

In addition to the parameter L, effective volumes may 
be further decreased depending on the value of the galactic 
convecting wind V c : the effective halo size L* = 2K(E)/V C 
plays a similar role as L (exponential cut-off of the con- 
tributions, [^lliEAMauri^ Hill). At low energy, L~L, 
decreasing V e g. This set the effective volume for p, 

^ ff = 2^ 2 {min(£,L*)} 3 . (A.I) 

Because of energy losses, the typical distance travelled 
by p ositr ons is 7% = \/AK (E)ti oss (e.g. IMaurin fc Tailletl 
120031: see lBrun et al.ll2007l and App. lB.2l for a more precise 
description). If r* 3> L, we recover V£ s , but if r* < L, the 
effective volume for positrons is 

V£ = \^frl (A.2) 

Note that the above volumes are in practise distorted in 
various directions (also because of the intrinsic spatial dis- 
tribution of DM sources), but they suffice for a qualitative 
estimate. 



E k 0.5 GeV 1 GeV 10 GeV 100 GeV 

A/ th (xl0 5 M Q ) p\e+ p\e + p\e+ p\e + 

max 340. /37. 340. /21. 340. /3.3 340. /0. 51 

med 0.32/1.6 1.2/1.2 6.4/0.42 6.4/0.15 

mm 0.0007/0.08 0.003/0.06 0.1/0.04 0.1/0.02 

Table A.l. Threshold mass (as defined in the text) in 
units of 10 5 M Q for p and e + at various energies for the three 
representative sets of propagation parameters (see Table [4] 
in Sect.SSl). 



Positrons: using r loss = 300 Myr x 1 GcV/E, and plugging 
the diffusion coefficient K(E) « K$E S in 7%, we get, at high 
energy, 

and 

For example, for the best propagation parameters set 
(Table ED, i-e. S = 0.7 and K = 0.0112, we get at 
E = 10 GeV, M t e h + ~ 5 x 10 4 M Q . 



A.2. Threshold mass Mth for the three propagation sets 

In Eqs. (|A.1[) and (|A.2j) , we set 77 = 10. We remind that this 
parameter sets the distance beyond which sources can be 
discarded. Taking a high value for 77 gives a conservative es- 
timate of M t h- The effective volumes, which do depend on 
energy, completely determine the number of clumps N c ff 
contributing to the flux. As fluctuations in N e f{ generate 
fluctuations in the signal, the mass threshold M t h is ob- 
tained demanding that N e g[M t h — 10M t h] > 10. We recall 
that for antimatter primaries, the most relevant contribu- 
tion comes from the local neighbourhood, so that the local 
number density of sub-halos dP(Ro)/dV can be used in 
the next approximations. 

Anti-protons: the effective number of clumps is given by 

dV 

N p eS [M - 10M] « JVtot [M - 10M] x — (R Q ) x V cS 
\M Q J \ lkpc / 

leading to 

<- 1 ° 4M - x {^rfef 1 } 3 - ( A - 3 ) 

For example, for the reference halo size L = 4 kpc and no 
galactic wind, we get the energy independent result MFL ~ 
6.4 x 10 5 M fl . 



" An analysis of these effective volumes has been prese nted 
and discussed in great details in IMaurin fc Tailletl (|2003l ), for 
both p and e + . 



Comparison: for the three sets of propagation parameters 
(as reminded in Table H] of Sect. 14. 3j) and for various en- 
ergies, Table |A~T1 gathers the threshold mass for p and e + , 
calculated from Eqs. (|A.3p and (|A.4|) . For anti-protons, we 
repeat that the leading parameter is L, except at low en- 
ergy when the Galactic wind blows particles efficiently out 
of the diffusive volume. If the halo size L is large {max 
set, i.e. L = 15 kpc), the variance associated with the flux 
of a clump mass range is sizable only for masses above 
> 3.10 7 M Q , independent of the energy. A small halo size 
(min set, i.e. L = 1 kpc) is associated with a strong wind, 
for which the threshold mass decreases down to ~ 10 2 M Q 
at 0.5 GeV (unmodulated). For positrons, the leading pa- 
rameter is the value of the diffusion coefficient, hence Kq 
and S. As for anti-protons, this is the min propagation set 
that leads to the smaller threshold mass. The range spanned 
by the various configurations is, however, tighter than that 
for anti-protons. The smaller values are observed at high 
energy (strong energy losses). At low energy and for the 
min set of parameters, as emphasised above, taking into 
account the wind would give extremely small M t h as for 
anti-protons. 

From these numbers, we may already predict that prop- 
agation para meters corresponding to small V e g will lead to 
smaller fluxes dBarrau et alJ2002Ll2005t iDonato et all2004t 
IMaurin et al.ll2006bh . but also a larger associated variance 
than the configurations with large V e g. 

A. 3. Choosing a higher M t h for the MC simulation 

In principle, in the MC realisations (Sect. 15. 3|) . all sub- 
halos with a mass larger than Mth need to be generated. 
We just provided an absolute criterion in App. IA.2I Using 
a relative criterion, a higher Mth (useful for reducing the 
computational time) can be found. 
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Mass (M so!ar ) 



Fig. A.l. Relative error cr/(</>) as a function of mass for 
several energies. The three lines are fitting results using 
function AM 13 , with (3 = 0.455, 0.445 and 0.45 from the 
top to bottom, respectively. The quantity \{)( 2l+1 )/ 2 is used 
to represent the mass bin Mj. 



The mean flux from clumps in the mass range [M, M + 
dM] is given by the number of clumps N(M) in the diffusive 
volume, and the variance of the signal is just proportional to 
V~N, so that the variance-flux-ratio is (see also lLavalle et al] 
120071) 



Using the relation N(M) cx il/ 1 ^" of Eq. dHJ), we get 



(A.5) 



££(M) cx M^— 1 )/ 2 

\4>) 



(A.6) 



thus recovering the results of lLavalle et all (|2007l ). who used 
sub-halos of equal masses and properties. As discussed in 
Sect. HOI for a m = 1.9, the luminosity is roughly constant 
per logarithmic mass bin, so that considering the sub-halo 
mass range Mj = [10 ? — lO t+1 ]M0, we get 



(70 (Mi) cx Al 



0.45 



(A.7) 



This equation states that the clumps maximally contribut- 
ing to the variance are the heaviest ones, and only these 
need to be taken into account. If we adopt M t h = 10 7 M Q , 
the neglected part will contribute only 0.1% to the total 
variance. 

As an illustration, Eq. (|A.6p is shown to be in full agree- 
ment with the result of the MC simulation in Fig. IA.1I (for 
anti-protons and a m = 1.9). The latter graphs are inde- 
pendent of the intrinsic profile of the clumps and of the 
propagation parameters. In Sect. 15751 we take advantage of 
this higher value to speed up the MC calculations. 

Appendix B: More details on mass- and 
space-related effects on fluxes and boosts 

We recall that the total contribution of sub-halos to the 
antimatter flux has a mean value given by: 



j tot 



4/T 



SxN cl x (Om x (G) v 



(B.l) 



which can be expressed as the product of two main terms. 
The first one is the integrated luminosity of clumps: 



(B.2) 



which is purely mass-related in the sense that it depends 
only on the mass distribution and internal features of sub- 
halos. The relevant parameters are Af m j n , a m , the mass- 
concentration relation and the inner sub-halo profile. The 
second is merely the propagation term {Q)v averaged over 
the spatial distribution of sub-halos, which is then of purely 
space type. 

B.l. Mass- related effects 

A simple analysis of the behaviour L c \ is helpful, and 
gives already interesting insights on the final results. In 
the limit for which M m ; n <C Af max and a m ~ 2, we find 
that the total number of clumps N c \ cx M^f" 1 AI m { n 1 ~ Qm 
(see Eq. [TH]) , and that the mean value (£)m cx M m i n Qm_1 x 
ln(M max /M min ) ~ Af min Qm_1 (assuming £ cx M c i), at vari- 
ance of small factors. This means that L c \ cx M"™ x 
ln(A/ max /A/ m i n ) increases with a m , and slightly (logarith- 
mically) depends on Af m j n . Actually, we find that spanning 
(a m ,M min ) respectively from (2,10^ 6 M Q ) to (1.8, 1O 6 M ) 
makes a decrease of only ~ 40 in the integrated luminosity 
(the above approximation gives ~ 150, but we recall that 
the actual dependence of £ in M c \ is not merely linear). 
Some numerical values are given in Table IB. II 

The pure mass-induced relative fluctuations of the 
sub- halo flux are given by ct£/(VN c i(£)m) (see Eq. |63|) . 
As c7£ cx Mjaiii m from the same arguments as 

above, those relative uncertainties approximately scale like 
(ln(M max /M min )) _1 , and are thus expected to only slightly 
(logarithmically) decrease when M m i n increases. They are 
actually found to lie in the range 10-0.1% for (a m ,M m i n ) 
going respectively from (2, 10~ 6 M Q ) to (1.8,10 6 M Q ) (see 



Table IBTT)) . This may appear surprising because we would 

1/9 

naively expect the relative variance to scale like N c \ ' cx 
M m i n ' am ~ 1 ^ 2 , and then to depend much more strongly on 
M m in- To summarise, the dropping of the total number of 
clumps, which reduces the statistical sample, is compen- 
sated by smaller fluctuations around the mean luminosity 
(£) from clump to clump (the range M m i n -M max gets thin- 
ner) , so that the mass-induced relative uncertainties remain 
roughly constant. 

B.2. Space-related effects for positrons 

The space-related effects for positrons are characterised by 
the mean value and the variance of Q over the spatial sub- 
halo distribution. For positrons, the relevant scale is the 
energy loss scale, w hich sets the charact eristic propagation 
length, as stressed in lLavalle et al.l (|2007f ). This propagation 
length is given by the following equation: 



Ads 



\ 



,5-1 



AKqt e 



1-5 



(B.3) 



where Kq and d are the normalisation and the logarith- 
mic slope of the diffusion coefficient, respectively, and e = 
(E/{Eq = 1 GeV}). If we take the medium propagation 
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Mmin 

(M e ) 


N c i 


(Qm 
(kpc 3 ) 




AT cl x (Om 
(kpc 3 ) 


lfr 6 


(1.6|4.1|1.0) x lO 13 ' 14 ! 16 


(85|10|3) x lO" 11 


(2.1|3.2|1.9) x 10 5 


(1.4|4.1|30.0) x 10 4 


1 


(2.6|1.6|1.0) x lO 8 ' 9 ! 10 


(49|16|7) x 1(T 6 


(9.1|10.1|7.7) x 10 2 


(1.3|2.5|7.4) x 10 4 


10 6 


(4.1|6.4|10.1) x 10 3 


(1.8 1.3|1.0) 


6.0|6.0|5.7 


(7.5 8.3 9.9) x 10 3 



Table B.l. Total number of clumps N c i, mean effective annihilation volume (£)m and its relative variance <7£/(£)m, and 
product N c \ x /(£) m (proportional to the total number of primary GCRs produced in clumps) for different mass models. 
We choose three different mass ranges varying the minimal mass from 10~ 6 to 10 6 A/q, and for different logarithmic 
slopes a m = {1.8|1.9|2.0} of the mass distribution. 



parameters of Table and a typical timescale for energy 
loss of te ~ 10 16 s, then for a 200 GeV injected energy, 
we find a propagation length A^ ~ 6.9 kpc x V e~ - 3 — 0.2, 
which ranges from 0.4 kpc at a detected energy Ed of 190 
GeV to 5.7 kpc at 1 GeV. This shows that for positrons, 
the main contributions to the flux arc likely to come from 
regions close to the solar neighbourhood. From Eq. (IB. 31) . 
we also see that a larger diffusion coefficient or/and a lower 
detected energy will allow the integration of contributions 
over a larger volume (the former case is, however, generally 
associated with a smaller diffusion slab model, which erases 
those extra contributions). 

Assuming now that all relevant contributions arc those 
inside a volume Vd around the Earth bounded by A^> < L 
(L is half the vertical extension of the diffusive halo), and 
that the propagation is roughly constant over this small 
volume, we can simplify the propagator Q e in the limit of 
infinite 3D diffusion: 



te 9(Xd 

X 



(7TA 2 D )3/2 



(B.4) 



Assuming also that the spatial distribution of clumps does 
not vary that much within Vd, thus given by dV(Ro)/dV , 
we get: 



(S 



, 3 ,+ dP 



X^(R Q ). (B.5) 



The interpretation is trivial, as well as the consequences for 
the boost estimate. Besides, the relative variance straight- 
forwardly reads: 



-1/2 



(B.6) 



This means that the relative variance decreases when the 
spatial probability function of clumps raises and when the 
effective detection volume increases, which is physically ob- 
vious but better quantified with the previous equations. If 
we argue in terms of local number density of clumps, this 
only says that the global space-related variance scales like 

— 1/2 

iV obs ' , where N Q ^ S is the number of clumps inside Vd- 



A more quantitative information is given in Table IB. 21 
In this table, we calculate the mean value and the vari- 
ance of Q e+ as defined by Eqs. fCTj) and respectively, 
for the three propagation models detailed in Table 0J We 
consider the injection of 200 GeV positrons in sources, 
(Q(Es) = 6(Es — 200 GeV)), and compute the propagator 
averaged on the spatial distribution of sub-halos together 
with the associated statistical variance. We show the results 
obtained for detected energies of 150 and 10 GeV, which 



correspond to diffusion lengths Ad of ~0.9 and 3.8 kpc, re- 
spectively. Such quantities are parts of the flux probability 
function related to a single clump, but encoding only the 
spatial and propagation information. 

We see from this table the expected behaviour when 
varying the propagation model: at a given detected energy, 
(Q e ) increases from the minimal to the maximal propa- 
gation configuration, and also increases when the positron 
is detected at a lower energy (its mean free path is some- 
how longer). There is a factor of ~ 20 between the min- 
imum (min model, high Ed) and maximum values (max 
model, low Ed). The space-associated contribution to the 
relative variance affecting the single clump flux is also given 
in Table [B~2l and is in the range 10-1000. It has to be com- 
pared with the relative mass-induced variance of Tabic [BTTT 
i.e. that on £. We see that while the relative variance on Q e+ 
is large, the one affecting £ almost always dominates, un- 
less the minimal mass of clumps is > 10 6 M o . Thus, though 
propagation uncertainties are important, the mass- induced 
effects are likely to outclass the systematic uncertainties 
over a large energy range. Nevertheless, they re-enter the 
game as soon as the propagation scale gets very short (de- 
tected energies very close to the injected energy). As the 
crossing space-mass term dominates the global relative vari- 
ance, we can determine the systematic errors affecting the 
global flux predictions, taking the previous ranges obtained 
for mass- like contributions. For the average clumpy contri- 
bution to the flux, we get ~ 40 x 20, which are three orders 
of magnitude. For the associated relative variance, we find 
ranges 0.1 — 10% (clump mass and number) and 10-1000 
(space-induced), which yield a total of four orders of mag- 
nitude. This provides the systematic uncertainties on the 
flux and its variance. Nevertheless, such uncertainties are 
diluted for the boost factor estimate, as only a small part of 
the parameter space gives a sub-halo contribution greater 
than that of the smooth. 



B.3. Space-related effects for anti-protons 

The same reasoning used for positrons can apply to anti- 
protons, that is the use of the propagation effective volume. 
Therefore, complementary to the following discussion, we 
refer the reader to the arguments and conclusion of the 
previous paragraph. 

The typical diffusion length for anti-protons depends 
mainly on the convective wind, and can be expressed as: 



Ad = 



_K(E) 



(B.7) 



where K(E) is the diffusion coefficient at energy E and 
V comr is the velocity of convection. This is quite different to 
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Propagation 


(s.kpc- 3 .GeV _1 ) 


is?) °gw/{G P ) 
(s.kpc- 3 .GeV _1 ) 


Mm 


5.462 x 10 b | 1.723 x 10 5 287.1|1089 


2.335 x 10 s |4.850 x 10' 829.2|622.2 


Med 


2.840 x 10' |1.744 x 10 b 67.95|446.6 


9.220 x 10 s 2.320 x 10 s 106.2|84.6 


Max 


3.666 x 10' |1.742 x 10 a 25.84 267.3 


2.652 x 10 y 1.157 x 10' J 18.42|15.22 


Med* 


5.589 x 10 s 3.191 x 10° 15.02 104.3 


1.801 x 10 1U |4.561 x 10 a 23.34 18.53 



Table B.2. Mean value and variance of Q (single clump) for positrons (left-hand side) and anti-protons (right hand side), 
each line accounting for each set of propagation parameters of Table [^J The results for positrons correspond to particles 
injected at 200 GeV in sources, and detected at 10| 150 GeV. For anti-protons, a flat spectrum dN/dT = 1 GcV 1 is 
injected at sources and detected kinetic energies of 10 1 150 GeV are also considered. For both positrons and anti-protons, 
the DM configuration is the reference one, except for Med*, for which we take a space distribution of sub-halos that tracks 
the smooth NFW component. 



that used for positrons in the sense that this length is much 
lower at low energy. With the medium set of propagation 
parameters, we get ~ 1.4|86 kpc for kinetic energies of 
0.1 1 10 GeV, respectively. However, the vertical boundary 
L of the diffusive halo limits that range to a few times 
L, so that the actual characteristic propagation length for 
anti-protons is usually comparable with the size of the slab. 
Anyhow, this means that above a few GeV, anti-protons can 
almost probe the entire diffusive slab, as they can originate 
from far away regions. Hence, the picture for anti-protons is 
the reversal of that for positrons. Besides, the characteristic 
diffusion length is larger for anti-protons than for positrons, 
so that the arguments using local quantities are less relevant 
here, unless for asymptotic values of the boost factor (which 
occur at low energy for anti-protons, see Eg. [68]) . 

Numerical values for the mean value and variance of Q p 
are given in Table IB. 21 where we have taken a flat injec- 
tion spectrum for anti-protons, dN/dT = 1 GeV -1 , and 
we have considered two detected kinetic energies of 10 and 
150 GeV (no energy losses for anti- protons). We recover the 
same range for systematic uncertainties as for positrons (see 

App.E2). 



